------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:/Users/silvio/Documents/CVR/ryp/output/s2_dokrigmi.log
  log type:  text
 opened on:   4 Jul 2018, 12:55:31

. scalar nstrata=58                      /*=58, CVR strata*/

. scalar nmiver=20

. 
. forval perpel= 1/3{
  2. scalar perpeh=`perpel'
  3. 
. matrix x2=J(nstrata, 1, 0)
  4. matrix modelmat_m=J(7, 1, 0)
  5. 
. 
. scalar diestim_mi=0
  6. scalar nestim_mi=0
  7. 
. local miper=1
  8.  while `miper'<=nmiver{
  9.  scalar miver=`miper'
 10. 
. use "${data}/${datafilemi}" , replace
 11. 
. gen j=i
 12. sort i perpe k
 13. 
. gen yy= y`miper'
 14. *ren y yy
. mkmat yy
 15. 
. 
. run "$code/cvrkrigrun"
 16. di miver
 17. matrix list sum_mao
 18. matrix list mao
 19. 
. matrix modelmat_m=modelmat_m+modelmat
 20. scalar diestim_mi=diestim_mi+diestim
 21. scalar nestim_mi=nestim_mi+nestim
 22. 
. 
. 
. 
. if miver==1{
 23. forval i = 1/58{
 24. matrix x2[`i',1]=mao[`i',2]*mao[`i',2]
 25. }
 26. matrix mao_mit=mao
 27. }
 28. else{
 29. matrix mao_mit=mao_mit+mao
 30. forval i = 1/58{
 31. matrix x2[`i',1]=x2[`i',1]+mao[`i',2]*mao[`i',2]
 32. }
 33. }
 34. 
. local miper=`miper'+1
 35. *drop yy
. } /*miper*/
 36.  *MI*********************************************************************************************************************
> *******************************
. 
. matrix mao_mit=mao_mit/nmiver 
 37. matrix x2=x2/nmiver 
 38. scalar madj=(nmiver+1)/(nmiver-1)
 39. 
. scalar diestim_mi=diestim_mi/nmiver 
 40. scalar nestim_mi=nestim_mi/nmiver 
 41. 
. 
. scalar jhat=1
 42. while jhat<=nstrata {
 43.         
.          matrix mao_mit[jhat,3]=mao_mit[jhat,3]+madj*(x2[jhat,1]-mao_mit[jhat,2]*mao_mit[jhat,2])            /*vnctotal*/   
>                      
 44.          matrix mao_mit[jhat,1]=int(mao_mit[jhat,1]+0.5)
 45.      matrix mao_mit[jhat,2]=int(mao_mit[jhat,2]+0.5) 
 46.          
.    scalar jhat =jhat+ 1
 47. } /* end loop */
 48. 
. 
. 
. *MI*************************************************************************************************************************
> ***************************
. 
. matrix ones58==J(1,58,1)
 49. matrix sum_mao_mit=ones58*mao_mit
 50. matrix list sum_mao_mit
 51. 
. di sqrt(sum_mao_mit[1,3])
 52. matrix list mao_mit
 53. di diestim
 54. di nestim
 55. 
. 
. if perpeh==1{
 56. matrix mao_mit_est=mao_mit
 57. matrix sum_mit_est=sum_mao_mit
 58. matrix modelmat_m_e=modelmat_m/nmiver 
 59. scalar diestim_mi_e=diestim_mi
 60. scalar nestim_mi_e=nestim_mi
 61. }
 62. else if perpeh==2{
 63. matrix mao_mit_sl=mao_mit
 64. matrix sum_mit_sl=sum_mao_mit
 65. matrix modelmat_m_s=modelmat_m/nmiver 
 66. scalar diestim_mi_s=diestim_mi
 67. scalar nestim_mi_s=nestim_mi
 68. }
 69. else if perpeh==3{
 70. matrix mao_mit_ot=mao_mit
 71. matrix sum_mit_ot=sum_mao_mit
 72. matrix modelmat_m_o=modelmat_m/nmiver 
 73. scalar diestim_mi_o=diestim_mi
 74. scalar nestim_mi_o=nestim_mi
 75. 
. }
 76. 
. 
. }
1

sum_mao[1,3]
         c1       c2       c3
r1    11122    27696  4401805

mao[58,3]
          c1       c2       c3
 r1      430     1087    20102
 r2       44       50       25
 r3       37       46       30
 r4       77      182     3339
 r5       89      271     6179
 r6       24       27        7
 r7       49      133     8736
 r8       26       29        5
 r9      213     1149   496080
r10       19       20        2
r11      798     2062    60296
r12      190      838    84053
r13       20       59     1126
r14     1608     3735    65384
r15      110      354     8614
r16       29      141    15120
r17       57      109      714
r18       21       31       45
r19      180      761   209125
r20      390     1061    48514
r21       76      119     1413
r22       70       73        6
r23      213      401     3137
r24       60       68       95
r25      973     1465     3655
r26      180      336     3282
r27       92      122      465
r28      162     1658  2374152
r29       52       86      270
r30      102      137      406
r31      179      598    29060
r32      403      691     6926
r33      101      218    11700
r34       66      143     1376
r35      493      839     8835
r36      648     1537    66513
r37      524      887     8986
r38       18       37      228
r39       33       34        2
r40      151      197      179
r41      109      307    25047
r42       63       76       41
r43       83      171     3667
r44       44       95     1211
r45      435      874    14144
r46       75      223     4707
r47      129      228      871
r48      170      256      524
r49       24       61      741
r50       39       65     1404
r51      154      201      402
r52      220      929   203657
r53      106      127       63
r54       23       32       39
r55       56       68       60
r56       73      110      537
r57      220     1650   465608
r58       92      432   130900
2

sum_mao[1,3]
         c1       c2       c3
r1    11149    27734  4349220

mao[58,3]
          c1       c2       c3
 r1      430     1087    20102
 r2       43       49       24
 r3       39       49       36
 r4       77      182     3339
 r5      100      319     8758
 r6       24       27        7
 r7       51      159    14040
 r8       27       30        6
 r9      211     1138   486675
r10       19       20        2
r11      806     2097    62823
r12      190      849    86838
r13       20       60     1200
r14     1614     3745    65610
r15      104      315     6641
r16       31      159    19584
r17       55      105      647
r18       21       31       45
r19      185      788   225522
r20      386     1034    45352
r21       73      108      980
r22       71       74        6
r23      213      401     3137
r24       62       71      103
r25      977     1474     3704
r26      188      353     3630
r27       97      129      521
r28      161     1635  2305336
r29       50       82      237
r30      103      139      415
r31      174      570    26136
r32      403      691     6926
r33      101      218    11700
r34       66      143     1376
r35      490      832     8639
r36      656     1563    69142
r37      524      887     8986
r38       17       33      181
r39       33       34        2
r40      148      192      167
r41      112      324    28500
r42       60       71       34
r43       85      181     4320
r44       44       95     1211
r45      426      832    12310
r46       75      223     4707
r47      120      204      654
r48      175      267      586
r49       24       61      741
r50       37       61     1200
r51      156      205      429
r52      223      948   212845
r53      104      124       56
r54       23       32       39
r55       56       68       60
r56       77      117      617
r57      219     1582   425256
r58       93      467   157080
3

sum_mao[1,3]
         c1       c2       c3
r1    11119    27630  4566272

mao[58,3]
          c1       c2       c3
 r1      431     1092    20345
 r2       43       49       24
 r3       38       48       33
 r4       82      204     4416
 r5      102      300     6600
 r6       24       27        7
 r7       52      172    17160
 r8       25       27        4
 r9      213     1158   505575
r10       19       20        2
r11      800     2063    60258
r12      194      880    93951
r13       20       61     1225
r14     1594     3677    62859
r15      104      315     6641
r16       30      150    17280
r17       55      105      647
r18       21       31       45
r19      185      788   225522
r20      394     1088    51783
r21       75      115     1260
r22       70       73        6
r23      214      403     3174
r24       60       68       95
r25      972     1461     3624
r26      182      340     3367
r27       90      120      443
r28      165     1727  2586672
r29       50       82      237
r30      102      137      406
r31      175      576    26708
r32      406      701     7232
r33      100      216    11406
r34       66      143     1376
r35      491      826     8379
r36      652     1530    65301
r37      519      867     8312
r38       18       37      228
r39       33       34        2
r40      149      193      171
r41      109      307    25047
r42       61       73       36
r43       84      176     3987
r44       44       95     1211
r45      435      867    13781
r46       73      214     4325
r47      131      234      930
r48      174      264      566
r49       25       65      878
r50       37       61     1200
r51      158      207      443
r52      220      890   183580
r53      102      120       49
r54       21       28       27
r55       54       65       52
r56       71      107      499
r57      211     1499   380365
r58       94      454   146520
4

sum_mao[1,3]
         c1       c2       c3
r1    11088    27507  4170965

mao[58,3]
          c1       c2       c3
 r1      436     1120    21656
 r2       43       49       24
 r3       38       48       33
 r4       82      208     4650
 r5       99      280     5608
 r6       24       27        7
 r7       48      128     7920
 r8       27       30        6
 r9      211     1147   496080
r10       19       20        2
r11      795     2081    63750
r12      187      805    76724
r13       24       85     2556
r14     1600     3693    63465
r15      106      330     7392
r16       31      159    19584
r17       53      193    25200
r18       20       29       37
r19      182      772   215610
r20      383     1013    43050
r21       75      115     1260
r22       70       73        6
r23      210      394     3029
r24       59       67       91
r25      972     1466     3700
r26      185      347     3497
r27       94      125      487
r28      159     1589  2170740
r29       50       82      237
r30      104      140      424
r31      175      576    26708
r32      402      684     6668
r33      102      221    11998
r34       64      138     1263
r35      497      848     9099
r36      649     1542    67165
r37      522      878     8651
r38       17       33      181
r39       33       34        2
r40      152      198      183
r41      110      313    26173
r42       59       70       32
r43       81      161     3067
r44       44       95     1211
r45      437      866    13670
r46       72      210     4140
r47      124      214      746
r48      173      263      564
r49       24       61      741
r50       37       61     1200
r51      157      207      457
r52      217      910   194673
r53      103      122       53
r54       22       30       33
r55       53       63       45
r56       72      108      517
r57      209     1495   379076
r58       95      491   175824
5

sum_mao[1,3]
         c1       c2       c3
r1    11118    27643  4271328

mao[58,3]
          c1       c2       c3
 r1      429     1081    19841
 r2       44       50       25
 r3       39       49       36
 r4       81      200     4189
 r5       95      272     5349
 r6       25       29        9
 r7       46      106     4620
 r8       27       30        6
 r9      211     1147   496080
r10       19       20        2
r11      799     2097    64881
r12      193      868    90889
r13       24       85     2556
r14     1602     3728    65729
r15      106      328     7298
r16       30      150    17280
r17       56      107      680
r18       20       29       37
r19      186      794   228876
r20      385     1027    44578
r21       73      108      980
r22       70       73        6
r23      212      399     3101
r24       64       73      113
r25      971     1463     3668
r26      184      344     3454
r27       93      124      476
r28      160     1612  2237532
r29       50       82      237
r30      101      136      397
r31      175      576    26708
r32      405      698     7130
r33      101      218    11700
r34       64      138     1263
r35      489      830     8574
r36      651     1602    78932
r37      516      852     7827
r38       17       33      181
r39       33       34        2
r40      151      197      179
r41      108      302    23946
r42       65       79       45
r43       83      171     3667
r44       44       95     1211
r45      428      845    12888
r46       75      223     4707
r47      132      236      947
r48      174      265      575
r49       24       61      741
r50       38       63     1300
r51      155      201      389
r52      229      986   231828
r53      105      125       60
r54       22       30       33
r55       54       64       45
r56       76      115      597
r57      218     1572   419508
r58       91      421   123420
6

sum_mao[1,3]
         c1       c2       c3
r1    11118    27611  4104958

mao[58,3]
          c1       c2       c3
 r1      433     1104    20859
 r2       44       50       25
 r3       38       48       33
 r4       81      200     4189
 r5      102      358    13118
 r6       25       29        9
 r7       46      106     4620
 r8       26       29        5
 r9      211     1129   477360
r10       19       20        2
r11      796     2176    76254
r12      191      847    86034
r13       21       65     1430
r14     1607     3742    65849
r15       98      277     4967
r16       31      159    19584
r17       55      105      647
r18       21       31       45
r19      183      777   218889
r20      384     1020    43811
r21       74      112     1116
r22       72       75        7
r23      214      403     3174
r24       66       76      122
r25      973     1459     3581
r26      190      357     3720
r27       91      121      454
r28      157     1543  2040192
r29       50       82      237
r30      102      137      406
r31      174      570    26136
r32      399      674     6380
r33       98      211    10828
r34       65      140     1319
r35      489      813     7938
r36      651     1547    67493
r37      517      857     7993
r38       17       33      181
r39       33       34        2
r40      149      193      171
r41      104      280    19788
r42       63       76       41
r43       84      176     3987
r44       44       95     1211
r45      434      849    12933
r46       74      219     4514
r47      132      236      947
r48      174      265      575
r49       25       65      878
r50       37       61     1200
r51      154      201      402
r52      219      923   200640
r53      105      125       60
r54       22       30       33
r55       55       66       52
r56       80      122      681
r57      222     1572   418500
r58       97      541   219336
7

sum_mao[1,3]
         c1       c2       c3
r1    11153    27883  4287829

mao[58,3]
          c1       c2       c3
 r1      428     1076    19602
 r2       45       52       27
 r3       37       46       30
 r4       79      195     3967
 r5       96      298     7493
 r6       24       27        7
 r7       51      159    14040
 r8       27       30        6
 r9      212     1139   486675
r10       19       20        2
r11      798     2143    70910
r12      190      841    84916
r13       22       72     1770
r14     1622     3787    67567
r15      116      395    11041
r16       31      159    19584
r17       52      120     7140
r18       22       33       53
r19      185      788   225522
r20      380      993    40808
r21       78      126     1746
r22       69       71        3
r23      215      405     3210
r24       60       68       95
r25      973     1459     3581
r26      185      347     3497
r27       96      128      510
r28      160     1612  2237532
r29       51       84      253
r30      101      136      397
r31      175      576    26708
r32      401      680     6571
r33      101      218    11700
r34       67      146     1435
r35      488      827     8510
r36      654     1612    80084
r37      526      897     9327
r38       18       37      228
r39       33       34        2
r40      149      193      171
r41      110      313    26173
r42       60       71       34
r43       82      166     3360
r44       44       95     1211
r45      433      853    13167
r46       74      219     4514
r47      129      228      871
r48      172      259      536
r49       24       61      741
r50       39       65     1404
r51      156      204      415
r52      225      961   219083
r53      107      129       67
r54       22       30       33
r55       54       65       49
r56       74      112      556
r57      222     1613   442735
r58       90      410   116160
8

sum_mao[1,3]
         c1       c2       c3
r1    11129    27741  4487980

mao[58,3]
          c1       c2       c3
 r1      434     1109    21130
 r2       43       49       24
 r3       38       48       33
 r4       74      169     2765
 r5       95      320    10286
 r6       24       27        7
 r7       49      133     8736
 r8       25       27        4
 r9      210     1137   486675
r10       19       20        2
r11      800     2063    60258
r12      189      826    81313
r13       20       61     1225
r14     1624     3810    69296
r15       96      267     4544
r16       31      159    19584
r17       56      107      680
r18       20       29       37
r19      184      783   222193
r20      390     1061    48514
r21       73      108      980
r22       70       73        6
r23      213      401     3137
r24       65       74      117
r25      974     1469     3686
r26      179      334     3239
r27       91      121      454
r28      163     1681  2443980
r29       49       80      221
r30       97      130      362
r31      177      587    27872
r32      403      691     6926
r33      103      223    12300
r34       65      140     1319
r35      497      848     9099
r36      656     1542    66569
r37      522      878     8651
r38       17       33      181
r39       33       34        2
r40      153      200      188
r41      106      291    21817
r42       59       70       32
r43       81      161     3067
r44       43       93     1145
r45      430      858    13483
r46       74      219     4514
r47      133      240      983
r48      173      263      564
r49       24       61      741
r50       36       59     1104
r51      158      207      443
r52      227      974   225410
r53      104      124       56
r54       22       30       33
r55       54       64       45
r56       74      112      556
r57      219     1620   448000
r58       91      443   139392
9

sum_mao[1,3]
         c1       c2       c3
r1    11143    27696  4216813

mao[58,3]
          c1       c2       c3
 r1      437     1126    21922
 r2       43       49       24
 r3       40       51       40
 r4       80      192     3763
 r5       94      321    10376
 r6       24       27        7
 r7       50      149    11880
 r8       27       30        6
 r9      211     1147   496080
r10       19       20        2
r11      797     2087    64113
r12      189      824    80885
r13       26       98     3486
r14     1608     3751    66315
r15      105      320     6859
r16       30      150    17280
r17       53      123     7560
r18       20       29       37
r19      184      783   222193
r20      388     1047    46920
r21       75      115     1260
r22       70       73        6
r23      212      399     3101
r24       67       77      127
r25      976     1482     3826
r26      186      349     3541
r27       90      120      443
r28      160     1612  2237532
r29       51       84      253
r30       99      133      379
r31      176      581    27287
r32      407      704     7336
r33      105      228    12915
r34       65      140     1319
r35      494      841     8901
r36      655     1559    68810
r37      517      854     7835
r38       17       33      181
r39       33       34        2
r40      150      195      175
r41      111      318    27324
r42       62       74       38
r43       84      176     3987
r44       43       93     1145
r45      434      863    13651
r46       72      210     4140
r47      126      221      805
r48      179      275      631
r49       25       65      878
r50       36       59     1104
r51      155      204      429
r52      220      890   183580
r53      107      129       67
r54       23       32       39
r55       54       64       45
r56       75      114      576
r57      215     1540   402497
r58       92      432   130900
10

sum_mao[1,3]
         c1       c2       c3
r1    11132    27971  4557903

mao[58,3]
          c1       c2       c3
 r1      432     1098    20611
 r2       43       49       24
 r3       37       46       30
 r4       76      178     3142
 r5       91      280     6615
 r6       24       27        7
 r7       49      137     9504
 r8       28       32        7
 r9      213     1149   496080
r10       19       20        2
r11      796     2047    59208
r12      186      927   128638
r13       22       72     1770
r14     1619     3784    67567
r15      110      354     8614
r16       30      150    17280
r17       56      132     8892
r18       20       29       37
r19      185      788   225522
r20      387     1040    46133
r21       73      108      980
r22       69       71        3
r23      216      407     3247
r24       61       70       99
r25      973     1459     3581
r26      185      347     3497
r27       90      120      443
r28      162     1658  2374152
r29       50       82      237
r30       99      133      379
r31      179      598    29060
r32      403      687     6766
r33      104      226    12606
r34       63      135     1209
r35      495      844     8966
r36      654     1556    68480
r37      515      848     7673
r38       17       33      181
r39       33       34        2
r40      149      193      171
r41      110      313    26173
r42       60       71       34
r43       82      166     3360
r44       43       93     1145
r45      437      880    14413
r46       74      219     4514
r47      128      226      847
r48      173      263      564
r49       25       65      878
r50       37       61     1200
r51      156      204      415
r52      218      917   197645
r53      104      124       56
r54       21       28       27
r55       55       67       60
r56       72      108      517
r57      230     1739   518320
r58       94      479   166320
11

sum_mao[1,3]
         c1       c2       c3
r1    11136    27830  4350194

mao[58,3]
          c1       c2       c3
 r1      432     1098    20611
 r2       43       49       24
 r3       37       46       30
 r4       79      191     3752
 r5       94      293     7296
 r6       25       29        9
 r7       49      139     9900
 r8       26       29        5
 r9      214     1159   505575
r10       19       20        2
r11      809     2195    75161
r12      191      848    86476
r13       21       65     1430
r14     1612     3750    66004
r15      106      328     7298
r16       30      150    17280
r17       53      100      583
r18       21       48      424
r19      183      777   218889
r20      386     1034    45352
r21       74      112     1116
r22       70       73        6
r23      213      401     3137
r24       61       70       99
r25      972     1472     3777
r26      181      338     3324
r27       93      124      476
r28      161     1635  2305336
r29       50       82      237
r30      100      135      388
r31      175      576    26708
r32      399      674     6380
r33      101      218    11700
r34       63      135     1209
r35      496      846     9033
r36      656     1563    69142
r37      517      853     7801
r38       17       33      181
r39       33       34        2
r40      147      190      163
r41      111      318    27324
r42       65       79       45
r43       84      176     3987
r44       44       95     1211
r45      440      882    14437
r46       74      219     4514
r47      130      231      896
r48      174      264      566
r49       25       65      878
r50       37       61     1200
r51      159      209      457
r52      216      904   191723
r53      105      154      664
r54       22       30       33
r55       55       67       60
r56       76      115      597
r57      216     1570   418966
r58       94      479   166320
12

sum_mao[1,3]
         c1       c2       c3
r1    11161    27952  4418384

mao[58,3]
          c1       c2       c3
 r1      430     1087    20100
 r2       43       49       24
 r3       38       48       33
 r4       78      187     3543
 r5       95      330    11102
 r6       25       29        9
 r7       48      128     7920
 r8       27       30        6
 r9      209     1136   486675
r10       19       20        2
r11      797     2087    64113
r12      192      852    87230
r13       23       78     2145
r14     1622     3813    69618
r15       99      288     5443
r16       31      159    19584
r17       56      107      680
r18       20       29       37
r19      185      788   225522
r20      387     1040    46133
r21       78      126     1746
r22       70       73        6
r23      216      407     3247
r24       65       77      195
r25      974     1474     3763
r26      186      349     3541
r27       91      121      454
r28      161     1635  2305336
r29       49       80      221
r30      100      135      388
r31      180      603    29663
r32      401      685     6726
r33      103      223    12300
r34       65      140     1319
r35      500      855     9300
r36      651     1547    67493
r37      524      890     9136
r38       17       33      181
r39       33       34        2
r40      152      198      183
r41      108      302    23946
r42       59       70       32
r43       85      181     4320
r44       44       95     1211
r45      436      863    13544
r46       72      210     4140
r47      128      226      854
r48      170      254      506
r49       25       65      878
r50       40       67     1512
r51      154      200      389
r52      225      961   219083
r53      106      127       63
r54       22       30       33
r55       54       64       45
r56       77      117      617
r57      223     1683   485012
r58       93      467   157080
13

sum_mao[1,3]
         c1       c2       c3
r1    11153    27793  4282487

mao[58,3]
          c1       c2       c3
 r1      437     1126    21925
 r2       43       49       24
 r3       38       48       33
 r4       79      191     3752
 r5       95      264     4979
 r6       24       27        7
 r7       46      106     4620
 r8       27       30        6
 r9      215     1178   524835
r10       19       20        2
r11      801     2138    70339
r12      187      801    75969
r13       23       78     2145
r14     1611     3754    66315
r15      103      308     6324
r16       29      141    15120
r17       55      105      647
r18       22       33       53
r19      188      805   235657
r20      391     1068    49321
r21       75      115     1260
r22       71       74        6
r23      214      403     3174
r24       66       76      122
r25      975     1472     3717
r26      185      347     3497
r27       87      116      411
r28      159     1589  2170740
r29       50       82      237
r30      101      136      397
r31      176      581    27287
r32      409      711     7545
r33      102      221    11998
r34       65      140     1319
r35      492      837     8769
r36      654     1556    68480
r37      523      882     8817
r38       17       33      181
r39       33       34        2
r40      151      197      179
r41      110      313    26173
r42       60       71       34
r43       82      166     3360
r44       45       98     1280
r45      436      855    13176
r46       74      219     4514
r47      129      228      871
r48      171      257      525
r49       24       61      741
r50       36       59     1104
r51      153      200      402
r52      223      948   212845
r53      106      127       63
r54       22       30       33
r55       52       61       38
r56       73      110      537
r57      227     1686   485680
r58       92      432   130900
14

sum_mao[1,3]
         c1       c2       c3
r1    11156    27898  4369743

mao[58,3]
          c1       c2       c3
 r1      436     1120    21656
 r2       43       49       24
 r3       37       46       30
 r4       82      204     4416
 r5       95      293     7267
 r6       24       27        7
 r7       47      117     6160
 r8       28       32        7
 r9      214     1168   515160
r10       19       20        2
r11      792     2107    67969
r12      189      837    84053
r13       25       91     3003
r14     1631     3804    68039
r15      111      364     9209
r16       30      150    17280
r17       54      102      614
r18       22       33       53
r19      185      788   225522
r20      382     1012    43050
r21       75      115     1260
r22       71       74        6
r23      209      392     2993
r24       61       70       99
r25      971     1463     3668
r26      177      330     3156
r27       95      127      498
r28      160     1612  2237532
r29       49       80      221
r30      104      140      424
r31      177      587    27872
r32      402      688     6826
r33       99      213    11115
r34       67      146     1435
r35      490      823     8317
r36      654     1558    68810
r37      519      863     8151
r38       17       33      181
r39       33       34        2
r40      151      197      179
r41      108      302    23946
r42       63       76       41
r43       84      176     3987
r44       43       93     1145
r45      429      854    13353
r46       75      223     4707
r47      131      234      921
r48      175      267      586
r49       26       70     1027
r50       37       61     1200
r51      163      213      457
r52      220      890   183580
r53      104      124       56
r54       22       30       33
r55       53       63       45
r56       76      115      597
r57      226     1696   492156
r58       94      502   185640
15

sum_mao[1,3]
         c1       c2       c3
r1    11173    28038  4439980

mao[58,3]
          c1       c2       c3
 r1      432     1098    20611
 r2       43       49       24
 r3       39       49       36
 r4       74      167     2650
 r5      100      292     6229
 r6       24       27        7
 r7       48      125     7392
 r8       28       32        7
 r9      213     1149   496080
r10       19       20        2
r11      799     2138    70434
r12      188      816    79310
r13       23       78     2145
r14     1632     3810    68355
r15      103      308     6324
r16       30      150    17280
r17       57      135     9360
r18       23       46      256
r19      186      794   228876
r20      388     1047    46920
r21       72      105      853
r22       70       73        6
r23      214      403     3174
r24       58       68      144
r25      975     1478     3795
r26      179      334     3239
r27       91      121      454
r28      160     1612  2237532
r29       50       82      237
r30      103      139      415
r31      177      587    27872
r32      405      693     6963
r33      102      221    11998
r34       66      143     1376
r35      489      830     8574
r36      652     1550    67821
r37      519      863     8151
r38       17       33      181
r39       33       34        2
r40      149      193      171
r41      108      302    23946
r42       63       76       41
r43       84      176     3987
r44       43       93     1145
r45      434      863    13651
r46       77      231     5106
r47      139      255     1134
r48      175      266      576
r49       24       61      741
r50       38       63     1300
r51      155      204      443
r52      221      936   206697
r53      104      124       56
r54       22       30       33
r55       56       68       60
r56       74      112      556
r57      231     1771   538692
r58       95      515   196560
16

sum_mao[1,3]
         c1       c2       c3
r1    11158    27880  4441030

mao[58,3]
          c1       c2       c3
 r1      435     1115    21389
 r2       44       50       25
 r3       39       49       36
 r4       80      192     3763
 r5       96      272     5330
 r6       24       27        7
 r7       48      128     7920
 r8       28       32        7
 r9      213     1158   505575
r10       19       20        2
r11      800     2072    61022
r12      191      848    86476
r13       21       65     1430
r14     1631     3809    68355
r15      113      374     9761
r16       31      159    19584
r17       54      126     7992
r18       20       29       37
r19      183      777   218889
r20      388     1047    46920
r21       71      101      735
r22       71       74        6
r23      215      405     3210
r24       65       74      117
r25      972     1466     3700
r26      184      344     3454
r27       91      121      454
r28      161     1635  2305336
r29       51       84      253
r30       99      133      379
r31      178      592    28463
r32      402      688     6826
r33       97      208    10545
r34       66      143     1376
r35      490      823     8317
r36      655     1562    69142
r37      519      863     8160
r38       18       37      228
r39       33       34        2
r40      151      197      179
r41      110      313    26173
r42       62       74       38
r43       82      166     3360
r44       44       95     1211
r45      427      835    12431
r46       74      219     4514
r47      126      220      792
r48      178      273      619
r49       25       65      878
r50       39       65     1404
r51      152      200      415
r52      215      865   172900
r53      105      125       60
r54       24       33       45
r55       56       67       56
r56       76      115      597
r57      221     1680   483555
r58       95      537   216580
17

sum_mao[1,3]
         c1       c2       c3
r1    11136    28559  7219729

mao[58,3]
          c1       c2       c3
 r1      436     1120    21656
 r2       43       49       24
 r3       39       49       36
 r4       84      213     4889
 r5       94      285     6782
 r6       24       27        7
 r7       47      119     6480
 r8       26       29        5
 r9      207     2025  3523284
r10       19       20        2
r11      809     2117    64339
r12      184      773    70091
r13       20       59     1126
r14     1610     3724    64613
r15      109      348     8341
r16       30      150    17280
r17       56      107      680
r18       19       27       30
r19      181      766   212355
r20      392     1075    50135
r21       73      108      980
r22       70       73        6
r23      213      401     3137
r24       60       68       95
r25      974     1463     3611
r26      184      344     3454
r27       97      129      521
r28      158     1566  2104960
r29       49       80      221
r30       96      129      354
r31      177      587    27872
r32      402      684     6668
r33       98      211    10828
r34       63      135     1209
r35      487      825     8446
r36      655     1559    68810
r37      517      857     7983
r38       17       33      181
r39       33       34        2
r40      155      203      196
r41      107      296    22869
r42       61       73       36
r43       83      171     3667
r44       44       95     1211
r45      433      867    13878
r46       74      219     4514
r47      129      229      879
r48      177      271      608
r49       25       65      878
r50       38       63     1300
r51      164      215      471
r52      216      866   172900
r53      103      122       53
r54       23       32       39
r55       53       62       37
r56       74      112      556
r57      230     1739   518320
r58       95      491   175824
18

sum_mao[1,3]
         c1       c2       c3
r1    11131    27889  4641882

mao[58,3]
          c1       c2       c3
 r1      427     1070    19352
 r2       44       50       25
 r3       39       49       36
 r4       79      191     3752
 r5       95      272     5349
 r6       25       29        9
 r7       49      137     9504
 r8       25       27        4
 r9      211     1129   477360
r10       19       20        2
r11      797     2118    68702
r12      188      825    81313
r13       19       55      990
r14     1615     3738    65217
r15      107      333     7537
r16       30      150    17280
r17       57      109      714
r18       21       31       45
r19      183      777   218889
r20      389     1066    49321
r21       75      115     1260
r22       71       74        6
r23      212      399     3101
r24       61       70       99
r25      973     1459     3581
r26      188      353     3630
r27       94      125      487
r28      163     1681  2443980
r29       50       82      237
r30      103      139      415
r31      177      587    27872
r32      402      684     6668
r33       99      213    11115
r34       64      138     1263
r35      492      828     8442
r36      649     1595    78168
r37      518      862     8151
r38       17       33      181
r39       33       34        2
r40      150      195      175
r41      109      307    25047
r42       61       73       36
r43       81      161     3067
r44       44       95     1211
r45      427      828    12099
r46       76      227     4904
r47      127      222      815
r48      175      267      586
r49       24       61      741
r50       37       61     1200
r51      161      213      485
r52      224      914   194580
r53      105      125       60
r54       23       32       39
r55       55       66       52
r56       73      110      537
r57      220     1668   476929
r58       99      617   295260
19

sum_mao[1,3]
         c1       c2       c3
r1    11152    28589  7442207

mao[58,3]
          c1       c2       c3
 r1      429     1082    19850
 r2       44       50       25
 r3       39       49       36
 r4       80      196     3967
 r5       91      276     6397
 r6       25       29        9
 r7       49      139     9900
 r8       27       30        6
 r9      210     2046  3593052
r10       19       20        2
r11      808     2107    63559
r12      188      825    81313
r13       25       91     3003
r14     1639     3810    67947
r15      103      310     6427
r16       29      141    15120
r17       54      126     7992
r18       21       31       45
r19      186      794   228876
r20      386     1034    45352
r21       73      108      980
r22       70       73        6
r23      212      399     3101
r24       60       71      158
r25      977     1480     3780
r26      184      344     3454
r27       91      121      454
r28      161     1635  2305336
r29       50       82      237
r30       96      129      354
r31      174      570    26136
r32      406      701     7232
r33      102      221    11998
r34       67      146     1435
r35      492      837     8769
r36      651     1549    67821
r37      518      858     7999
r38       17       33      181
r39       33       34        2
r40      149      193      171
r41      107      296    22869
r42       64       77       43
r43       81      161     3067
r44       43       93     1145
r45      431      847    12919
r46       74      219     4514
r47      130      231      896
r48      178      273      619
r49       24       61      741
r50       36       59     1104
r51      159      209      457
r52      224      873   173987
r53      103      122       53
r54       21       28       27
r55       54       64       45
r56       74      112      556
r57      220     1592   431043
r58       94      502   185640
20

sum_mao[1,3]
         c1       c2       c3
r1    11164    27932  4389808

mao[58,3]
          c1       c2       c3
 r1      438     1132    22195
 r2       43       49       24
 r3       38       48       33
 r4       81      200     4189
 r5       98      301     7646
 r6       24       27        7
 r7       49      137     9504
 r8       26       29        5
 r9      211     1129   477360
r10       19       20        2
r11      794     2122    69225
r12      190      836    83657
r13       21       66     1485
r14     1633     3814    68514
r15      111      362     9074
r16       29      141    15120
r17       55      105      647
r18       20       43      337
r19      186      794   228876
r20      382     1003    41918
r21       79      129     1925
r22       71       74        6
r23      212      399     3101
r24       63       72      108
r25      975     1472     3717
r26      178      332     3197
r27       95      127      498
r28      161     1635  2305336
r29       51       84      253
r30       99      133      379
r31      180      603    29663
r32      403      691     6926
r33       99      213    11115
r34       65      140     1319
r35      488      819     8193
r36      655     1559    68810
r37      521      880     8817
r38       18       37      228
r39       33       34        2
r40      149      193      171
r41      109      307    25047
r42       60       71       34
r43       82      166     3360
r44       45       98     1280
r45      430      851    13140
r46       75      223     4707
r47      124      214      746
r48      173      263      564
r49       24       61      741
r50       37       61     1200
r51      162      212      457
r52      222      902   189040
r53      106      127       63
r54       22       30       33
r55       56       69       68
r56       75      114      576
r57      227     1747   524270
r58       92      432   130900

sum_mao_mit[1,3]
           c1         c2         c3
r1      11143      27872  4772788.2
2184.6712

mao_mit[58,3]
            c1         c2         c3
 r1        433       1101  21140.642
 r2         43         49  25.055132
 r3         38         48  35.639474
 r4         79        192  3975.2205
 r5         96        295  8056.4857
 r6         24         28  8.5284211
 r7         49        133  9360.8462
 r8         27         30  8.2031579
 r9        212       1236  880129.14
r10         19         20          2
r11        800       2106  67457.941
r12        189        838  86558.672
r13         22         72  2030.2558
r14       1617       3764  68517.352
r15        106        329  8537.6232
r16         30        151  17701.625
r17         55        116  4593.7547
r18         21         33  123.37842
r19        184        784  223184.06
r20        387       1041  46845.305
r21         75        113    1260.36
r22         70         73  6.7005263
r23        213        401    3156.01
r24         62         71  126.07842
r25        974       1468  3742.4058
r26        184        343   3492.463
r27         92        123  483.40211
r28        161       1628  2288280.6
r29         50         82  241.25842
r30        101        135     405.59
r31        177        584  27708.308
r32        403        690  6975.9921
r33        101        218  11707.371
r34         65        141  1338.6153
r35        492        834  8784.2578
r36        653       1559   70112.56
r37        520        869  8597.2163
r38         17         34  196.06579
r39         33         34          2
r40        150        195  186.40934
r41        109        306  24978.319
r42         62         74  46.576184
r43         83        170  3679.2836
r44         44         95  1200.3216
r45        433        857  13582.897
r46         74        219  4576.3495
r47        129        228  992.74355
r48        174        265  604.92303
r49         25         63  817.00513
r50         37         62  1247.4821
r51        157        206  456.90053
r52        221        919  199846.45
r53        105        126  139.31105
r54         22         30  36.330395
r55         54         65  56.130921
r56         75        113  582.68566
r57        221       1636  464868.72
r58         94        477  169980.39
58
11164
1

sum_mao[1,3]
         c1       c2       c3
r1    10392    20547  1612172

mao[58,3]
         c1      c2      c3
 r1     193     384   18640
 r2       2       4       6
 r3      29      60    1124
 r4      62     126    4467
 r5     141     285   42660
 r6      19      39     445
 r7      18      38     491
 r8      10      21     152
 r9      58     120    4928
r10      21      40     395
r11     665     808   12691
r12      58     117    5658
r13      22      42     592
r14    1549    3070  333534
r15     149     298    7905
r16      57      68    2452
r17      60      96     992
r18      56      75    2060
r19     183     334   16531
r20     386     786   49095
r21      35      79     222
r22      65     129    1016
r23     185     356   12712
r24     152     297    9731
r25     433     997    9369
r26     242     542   16778
r27     149     332    4869
r28      89     183    2904
r29      99     260    1871
r30     160     340    8200
r31     267     589   20896
r32     182     846  278880
r33     204     547   13071
r34     117     116    4170
r35     287     343    3894
r36     464     694   39934
r37     197     479    5983
r38      13      35      34
r39      21      56      62
r40     237     692    9923
r41      33     111     143
r42      92     203     933
r43      34      89     176
r44       7      21       6
r45     452     810  105588
r46     129     295    1665
r47     298     571   11168
r48     343     569   42991
r49      10      28      17
r50      84     160    4538
r51     349     692   79119
r52     207     400   24385
r53      34      63     707
r54      10      18      66
r55      32      60     640
r56      83     154    3604
r57     842    1546  386642
r58      17      34     447
2

sum_mao[1,3]
        c1      c2      c3
r1   10360   17276  976103

mao[58,3]
         c1      c2      c3
 r1     193     249   13344
 r2       3       4      10
 r3      26      36     669
 r4      63      86    3514
 r5     130     243   27226
 r6      19      26     326
 r7      16      23     287
 r8       8      11      68
 r9      60      85    3626
r10      21      28     193
r11     652     792   12187
r12      57     113    6272
r13      21      35     207
r14    1543    1733   40666
r15     155     239    6343
r16      55      67    1722
r17      61      95     768
r18      56      74    1556
r19     177     264   10335
r20     388     695   36011
r21      37      75     185
r22      64     109     709
r23     185     303    8937
r24     150     259    6894
r25     429     991    9445
r26     235     494   11785
r27     144     310    3398
r28      90     168    2192
r29     101     258    1458
r30     159     316    6021
r31     272     547   16057
r32     184     856  285600
r33     203     509    9393
r34     117     120    3120
r35     290     347    3977
r36     457     683   38732
r37     198     438    4523
r38      14      35      29
r39      21      52      47
r40     242     682    7745
r41      30      98      88
r42      96     187     760
r43      31      75     110
r44       7      20       5
r45     460     574   16402
r46     129     260    1238
r47     307     608   13114
r48     338     450   16775
r49      10      26      13
r50      86     119    1676
r51     350     580   44467
r52     205     299   15514
r53      36      50     446
r54       9      12      26
r55      31      43     337
r56      79     116    2148
r57     845    1281  277361
r58      15      28      46
3

sum_mao[1,3]
        c1      c2      c3
r1   10380   16549  971837

mao[58,3]
         c1      c2      c3
 r1     194     240   12475
 r2       3       6       8
 r3      28      47     643
 r4      58      86    2579
 r5     129     208   27563
 r6      19      32     267
 r7      15      26     209
 r8      10      18      85
 r9      56     103    2524
r10      21      35     144
r11     661     803   12535
r12      52      73    3248
r13      21      27     373
r14    1562    1754   41580
r15     153     289    5149
r16      56      93    1576
r17      62     128     699
r18      56     102    1372
r19     178     228    8735
r20     380     468   25861
r21      36      41      61
r22      65      75     496
r23     184     221    7009
r24     151     196    5707
r25     434    1001    9475
r26     239     378    9476
r27     151     274    3099
r28      85     118    1568
r29     101     208    1106
r30     160     244    4921
r31     271     368   11255
r32     179     831  268950
r33     205     590    8431
r34     117     166    2752
r35     289     478   30188
r36     459     686   39073
r37     201     229     900
r38      13      22      14
r39      21      36      17
r40     241     525    4951
r41      33      93      74
r42      95     133     306
r43      33      50      29
r44       7      19       4
r45     449     560   15624
r46     131     284    1113
r47     295     546    9857
r48     339     562   41987
r49       9      16       4
r50      86     139    1461
r51     348     690   78660
r52     208     411   27834
r53      38      62     283
r54      12      18      30
r55      33      55     207
r56      85     136    1848
r57     849    1303  235228
r58      14      19     214
4

sum_mao[1,3]
         c1       c2       c3
r1    10431    19151  1918779

mao[58,3]
         c1      c2      c3
 r1     190     226   32763
 r2       3      11       3
 r3      28      78    1249
 r4      58     122    6324
 r5     132     263   79500
 r6      18      49     475
 r7      19      59     536
 r8       9      34      17
 r9      59     226   37968
r10      21      52     356
r11     669     813   12769
r12      59      93   11475
r13      18      25     755
r14    1567    1759   41688
r15     153     343   14601
r16      55     189    3463
r17      64     303   72720
r18      57     230    2765
r19     181     262   26154
r20     389     626   87414
r21      36      68     423
r22      65     101    1769
r23     188     279   22293
r24     153     237   17315
r25     435     993    9091
r26     238     456   29228
r27     146     292    8461
r28      91     152    5412
r29     101     244    3530
r30     158     283   14372
r31     271     492   38510
r32     184     852  282230
r33     202     768   19395
r34     119     258    6970
r35     282     337    3759
r36     461     803   68743
r37     199     416   11076
r38      14      33      71
r39      21      50     113
r40     238     639   18147
r41      32     103     244
r42      96     179    1845
r43      36      83     359
r44       7      21      11
r45     449     671   37352
r46     132     287    3131
r47     304     602   12841
r48     341     454   17075
r49      10      25      31
r50      86     189    3792
r51     351     810  123789
r52     210     478   35732
r53      34      66     927
r54      11      21      88
r55      34      67     933
r56      85     151    6027
r57     850    1410  680267
r58      12      18     432
5

sum_mao[1,3]
         c1       c2       c3
r1    10384    18203  1270033

mao[58,3]
         c1      c2      c3
 r1     194     356   14685
 r2       2       4       4
 r3      27      54     750
 r4      60     111    3314
 r5     136     228   35998
 r6      18      37     304
 r7      20      40     467
 r8       9      17      89
 r9      58     112    3531
r10      21      35     200
r11     666     809   12652
r12      55      94    4276
r13      18      31     323
r14    1551    2568  174190
r15     152     410    6333
r16      56      68    1852
r17      61      98     794
r18      57      76    1670
r19     177     262   10362
r20     390     511   32070
r21      38      43      80
r22      65      74     583
r23     185     239    8363
r24     147     215    6376
r25     435    1006    9581
r26     239     401   11137
r27     147     278    3450
r28      91     139    2114
r29     101     212    1300
r30     160     263    5785
r31     270     387   13135
r32     180     836  272240
r33     204     524    9808
r34     119     122    3344
r35     291     348    4004
r36     461     643   23515
r37     205     233     937
r38      13      23      16
r39      21      35      20
r40     239     529    5720
r41      34      95      92
r42      91     117     330
r43      33      49      34
r44       7      18       4
r45     457     570   16188
r46     129     246    1268
r47     294     544    9786
r48     340     452   16975
r49      10      18       6
r50      85     127    1693
r51     352    1038  237356
r52     199     450   15190
r53      34      64     412
r54      11      19      40
r55      34      64     420
r56      81     161    2329
r57     841    1678  282321
r58      13      22     217
6

sum_mao[1,3]
         c1       c2       c3
r1    10390    18723  1275472

mao[58,3]
         c1      c2      c3
 r1     186     329   12627
 r2       2       4       4
 r3      28      59     755
 r4      60     115    3100
 r5     128     228   29828
 r6      18      38     285
 r7      20      43     437
 r8      10      22     103
 r9      58     127    3303
r10      21      43     187
r11     670     833    6831
r12      56      98    4147
r13      21      36     411
r14    1551    2568  174190
r15     159     392    6482
r16      55      90    1671
r17      62     127     768
r18      56     101    1508
r19     180     274   10024
r20     391     514   30153
r21      37      42      71
r22      63      73     512
r23     184     243    7738
r24     146     210    5883
r25     433    1010    9883
r26     232     384    9816
r27     149     279    3316
r28      94     141    2110
r29     101     211    1216
r30     161     260    5479
r31     271     385   12378
r32     184     856  285600
r33     207     585    9446
r34     118     167    3076
r35     287     475   29767
r36     461     689   39416
r37     204     232     928
r38      14      24      17
r39      21      35      19
r40     241     528    5440
r41      38     106     108
r42      93     119     322
r43      33      49      32
r44       7      18       4
r45     452     787   65981
r46     130     254    1205
r47     295     565   10930
r48     340     452   16975
r49       9      16       5
r50      86     167    1621
r51     352     925  176841
r52     208     468   15523
r53      33      66     363
r54      11      21      37
r55      33      66     370
r56      77     152    1969
r57     841    1601  264088
r58      12      21     173
7

sum_mao[1,3]
         c1       c2       c3
r1    10347    19873  1236359

mao[58,3]
         c1      c2      c3
 r1     197     370   14629
 r2       1       2       1
 r3      29      54     836
 r4      62     115    3419
 r5     132     247   32759
 r6      19      35     328
 r7      15      27     254
 r8       9      16      86
 r9      59     102    3532
r10      21      32     194
r11     666     846    8530
r12      53      98    3836
r13      20      37     385
r14    1537    2745  221064
r15     142     310    5340
r16      55      62    1729
r17      64     100     847
r18      55      69    1506
r19     178     313   10475
r20     395     796   37404
r21      33      74     147
r22      66     128     756
r23     183     343    8764
r24     152     297    7095
r25     434    1015    9994
r26     237     535   12013
r27     144     325    3405
r28      90     186    2197
r29     100     263    1432
r30     161     346    6187
r31     271     598   15974
r32     184     852  282230
r33     204     527    9506
r34     116     113    3075
r35     291     348    4004
r36     459     640   23310
r37     195     474    4397
r38      13      35      25
r39      21      56      47
r40     239     698    7571
r41      32     108     101
r42      96     211     762
r43      35      91     140
r44       7      21       5
r45     452     564   15834
r46     130     294    1260
r47     295     546    9857
r48     341     566   42488
r49      10      28      13
r50      84     121    1603
r51     349     805  122360
r52     203     386   15271
r53      32      54     354
r54      11      17      38
r55      33      56     383
r56      81     149    2264
r57     839    1599  274063
r58      15      28     280
8

sum_mao[1,3]
         c1       c2       c3
r1    10363    17445  1271830

mao[58,3]
         c1      c2      c3
 r1     188     231   18264
 r2       3       5      14
 r3      27      41    1006
 r4      65      89    5070
 r5     134     232   39788
 r6      19      29     459
 r7      18      29     508
 r8      10      17     157
 r9      61     101    5630
r10      21      33     408
r11     663     805   12574
r12      57      86    5643
r13      21      26     557
r14    1535    1703   31380
r15     159     250    9296
r16      54      64    2273
r17      60      95    1024
r18      56      73    2127
r19     179     279   16333
r20     386     702   50699
r21      38      77     270
r22      65     112    1049
r23     184     309   12986
r24     146     256    9272
r25     432     992    9264
r26     243     513   17470
r27     150     323    5095
r28      88     166    2932
r29     102     260    2051
r30     165     329    9005
r31     269     545   21903
r32     184     856  285600
r33     202     515   13234
r34     118     119    4380
r35     282     337    3759
r36     455     680   38392
r37     199     441    6304
r38      14      35      40
r39      21      52      64
r40     236     664   10160
r41      36     118     175
r42      97     189    1071
r43      36      87     204
r44       8      23       8
r45     455     664  110490
r46     130     263    1746
r47     293     542    9715
r48     338     450   16775
r49      10      26      18
r50      87     138    5027
r51     346     573   43446
r52     201     314   23743
r53      36      55     819
r54      11      16      82
r55      34      52     746
r56      82     127    3633
r57     840    1316  397379
r58      14      21     313
9

sum_mao[1,3]
         c1       c2       c3
r1    10353    18653  1489809

mao[58,3]
         c1      c2      c3
 r1     186     399   14185
 r2       3       9       1
 r3      26      69     483
 r4      60     148    3055
 r5     133     264   36469
 r6      19      50     236
 r7      17      47     193
 r8       8      23       6
 r9      58     169   18928
r10      21      39     147
r11     665     808   12613
r12      56     113    4671
r13      16      33     270
r14    1544    3274  404675
r15     152     387    6257
r16      56      70    1961
r17      64     104     924
r18      56      78    1705
r19     179     238   10913
r20     387     485   33261
r21      36      41      76
r22      65      72     616
r23     186     225    8871
r24     145     197    6514
r25     430     982    9056
r26     237     388   11575
r27     151     280    3853
r28      91     132    2228
r29     100     211    1350
r30     164     260    6418
r31     266     372   13464
r32     179     831  268950
r33     200     521    9920
r34     118     124    3484
r35     286     342    3867
r36     456     682   38561
r37     204     237    1296
r38      14      25      20
r39      21      36      21
r40     238     530    6009
r41      31      87      81
r42      94     122     373
r43      33      50      37
r44       8      21       6
r45     453     565   15905
r46     131     252    1384
r47     299     554   10143
r48     335     446   16477
r49       9      17       6
r50      87     139    1770
r51     354     931  178895
r52     208     411   27834
r53      33      56     265
r54      10      16      26
r55      34      58     273
r56      82     142    2126
r57     845    1462  286840
r58      14      29     266
10

sum_mao[1,3]
         c1       c2       c3
r1    10374    18359  1473774

mao[58,3]
         c1      c2      c3
 r1     190     409   13700
 r2       3       6      11
 r3      29      61     869
 r4      65     141    3814
 r5     139     285   36518
 r6      19      39     341
 r7      18      37     384
 r8       8      16      76
 r9      57     111    3696
r10      21      36     300
r11     665     808   12691
r12      60     127    4948
r13      19      41     349
r14    1546    3386  444884
r15     148     298    5860
r16      56      71    1799
r17      61     100     776
r18      56      78    1570
r19     178     234   11365
r20     388     490   31776
r21      38      44      78
r22      66      74     610
r23     181     222    8279
r24     150     202    6612
r25     433    1010    9883
r26     236     388   10573
r27     151     283    3536
r28      88     127    1941
r29     101     214    1261
r30     162     257    5782
r31     266     374   12497
r32     182     846  278880
r33     200     523    9486
r34     120     128    3300
r35     285     341    3840
r36     460     688   39245
r37     204     237    1296
r38      14      25      18
r39      21      36      19
r40     241     539    5647
r41      32      90      80
r42      96     124     356
r43      35      53      38
r44       8      21       5
r45     447     665   77872
r46     130     251    1259
r47     298     571   11168
r48     340     452   16975
r49       9      17       5
r50      86     145    3607
r51     349     692   79119
r52     209     376   18676
r53      35      56     562
r54      10      16      50
r55      33      54     510
r56      85     134    2820
r57     834    1282  281951
r58      13      28     211
11

sum_mao[1,3]
        c1      c2      c3
r1   10362   16488  996626

mao[58,3]
         c1      c2      c3
 r1     192     305   12394
 r2       3       5       9
 r3      27      42     646
 r4      62     102    3049
 r5     134     245   30113
 r6      18      28     262
 r7      17      27     291
 r8      10      15      95
 r9      57      87    2938
r10      21      28     173
r11     656     797   12264
r12      55      95    3685
r13      21      36     378
r14    1547    2498  148394
r15     151     207    5385
r16      56      71    1596
r17      62     100     707
r18      55      76    1340
r19     180     212    9234
r20     390     523   27634
r21      35      49      58
r22      65      81     502
r23     185     230    7206
r24     151     202    5797
r25     435     980    8632
r26     242     415    9839
r27     148     280    3014
r28      90     133    1782
r29     101     224    1120
r30     162     264    5110
r31     268     409   11151
r32     187     867  292400
r33     204     524    8451
r34     120     127    2930
r35     283     339    3785
r36     456     682   38561
r37     203     300    4850
r38      14      28      16
r39      21      41      17
r40     243     582    5095
r41      31      91      66
r42      91     136     284
r43      33      59      30
r44       7      19       4
r45     447     558   15484
r46     129     252    1093
r47     294     563   10851
r48     339     451   16875
r49       9      19       5
r50      86     125    1493
r51     347     461   17582
r52     206     265   14026
r53      35      43     377
r54      10      12      28
r55      33      40     341
r56      80      96    1958
r57     844    1018  245009
r58      14      24     217
12

sum_mao[1,3]
         c1       c2       c3
r1    10322    19479  1527235

mao[58,3]
         c1      c2      c3
 r1     191     246   22543
 r2       3      11       2
 r3      28      82     851
 r4      61     133    4764
 r5     135     219   56621
 r6      18      52     323
 r7      17      55     292
 r8       9      34      11
 r9      59     226   37968
r10      21      51     243
r11     666     846    8530
r12      51      66    5839
r13      18      23     514
r14    1532    1699   31173
r15     154     428   10072
r16      55      65    2856
r17      61      97    1268
r18      57      76    2668
r19     177     330   17089
r20     389     743   59787
r21      33      68     243
r22      65     116    1209
r23     182     330   14286
r24     146     280   10781
r25     433     984    8889
r26     237     513   19814
r27     150     328    6099
r28      90     180    3620
r29     102     260    2460
r30     162     336   10327
r31     266     555   25377
r32     185     861  288990
r33     202     524   15301
r34     118     117    5258
r35     278     333    3652
r36     461     689   39416
r37     193     429    7112
r38      14      35      49
r39      21      52      77
r40     238     670   12398
r41      34     111     188
r42      95     185    1232
r43      32      77     194
r44       7      20       8
r45     450     673   37520
r46     131     266    2113
r47     299     554   10143
r48     339     451   16875
r49       9      24      17
r50      83     152    2436
r51     352    1038  237356
r52     204     553   22964
r53      33      76     595
r54      11      24      60
r55      34      79     636
r56      79     180    3550
r57     839    1857  452231
r58      13      17     345
13

sum_mao[1,3]
         c1       c2       c3
r1    10349    20102  1271384

mao[58,3]
         c1      c2      c3
 r1     184     343   12370
 r2       3       6       9
 r3      28      55     756
 r4      61     119    3208
 r5     134     259   32724
 r6      19      37     318
 r7      19      38     395
 r8       9      18      84
 r9      56     109    3085
r10      21      39     188
r11     661     802   12341
r12      58     109    4453
r13      18      34     302
r14    1546    2849  265118
r15     152     298    5930
r16      56      66    1738
r17      62      99     771
r18      55      73    1460
r19     175     318    9815
r20     384     780   34265
r21      36      81     170
r22      64     127     689
r23     184     353    8589
r24     146     284    6345
r25     432     992    9264
r26     238     532   11743
r27     154     342    3775
r28      91     187    2177
r29     101     265    1416
r30     160     339    5923
r31     270     594   15370
r32     178     826  265680
r33     203     540    9124
r34     118     117    3084
r35     287     343    3894
r36     459     686   39073
r37     198     481    4394
r38      14      37      28
r39      21      56      45
r40     239     698    7338
r41      32     108      97
r42      96     211     738
r43      35      91     136
r44       6      18       3
r45     450     783   65394
r46     130     297    1221
r47     298     571   11168
r48     341     566   42488
r49      10      28      12
r50      87     160    1667
r51     351     696   80040
r52     204     386   14949
r53      34      62     387
r54      11      20      37
r55      36      66     442
r56      83     152    2304
r57     834    1524  262502
r58      17      32     348
14

sum_mao[1,3]
         c1       c2       c3
r1    10362    17615  1324636

mao[58,3]
         c1      c2      c3
 r1     189     334   14784
 r2       3       9       1
 r3      29      73     607
 r4      57     133    2785
 r5     134     274   37377
 r6      19      47     238
 r7      19      51     244
 r8       7      20       5
 r9      55     160   16960
r10      21      45     192
r11     671     843    7678
r12      57     107    4885
r13      17      33     307
r14    1540    2945  295530
r15     148     237    5995
r16      56      75    1980
r17      63     105     907
r18      55      81    1663
r19     178     220   11770
r20     392     477   35315
r21      36      41      77
r22      64      71     625
r23     189     222    9716
r24     151     190    7278
r25     435     993    9091
r26     244     390   12420
r27     146     267    3638
r28      90     123    2221
r29     102     218    1417
r30     158     241    6037
r31     269     368   14015
r32     182     846  278880
r33     206     550   10886
r34     116     130    3398
r35     288     345    3921
r36     457     796   67546
r37     201     233    1257
r38      14      25      20
r39      21      36      21
r40     238     536    6067
r41      34      96      99
r42      92     119     361
r43      33      50      37
r44       8      21       6
r45     456     756   74259
r46     129     250    1362
r47     295     565   10930
r48     339     451   16875
r49       8      15       4
r50      85     165    3155
r51     343     456   17176
r52     209     413   28107
r53      36      62     365
r54      11      20      41
r55      35      62     334
r56      81     118    2205
r57     838    1081  291335
r58      13      25     231
15

sum_mao[1,3]
         c1       c2       c3
r1    10338    17857  1232023

mao[58,3]
         c1      c2      c3
 r1     186     324   12949
 r2       3       6      11
 r3      27      53     743
 r4      66     125    3878
 r5     131     242   31992
 r6      19      37     336
 r7      19      38     422
 r8       8      16      75
 r9      56     113    3519
r10      21      40     296
r11     666     808   12574
r12      59     105    4719
r13      19      34     344
r14    1536    2643  195137
r15     152     297    6097
r16      56      73    1774
r17      60     101     740
r18      54      78    1440
r19     177     252   11084
r20     386     511   31019
r21      39      46      81
r22      65      80     583
r23     183     241    8347
r24     152     212    6697
r25     431     987    9160
r26     243     398   11057
r27     150     277    3442
r28      90     132    2003
r29     101     211    1244
r30     159     252    5494
r31     269     383   12606
r32     181     841  275550
r33     203     561    9640
r34     116     123    3042
r35     290     347    3977
r36     461     689   39416
r37     203     236    1283
r38      14      25      18
r39      21      37      19
r40     241     534    5570
r41      34      96      89
r42      92     130     323
r43      33      50      33
r44       8      21       5
r45     448     777   77151
r46     127     276    1185
r47     286     547   10230
r48     339     562   41987
r49      10      18       6
r50      85     165    3476
r51     352     698   80503
r52     207     385   18069
r53      36      62     586
r54      11      19      59
r55      31      54     444
r56      81     135    2526
r57     832    1331  276765
r58      13      23     208
16

sum_mao[1,3]
         c1       c2       c3
r1    10338    20556  1505914

mao[58,3]
         c1      c2      c3
 r1     187     387   14521
 r2       2       4       5
 r3      27      56     798
 r4      60     122    3527
 r5     134     251   37190
 r6      19      40     361
 r7      18      37     403
 r8       8      15      75
 r9      57     109    3632
r10      21      35     213
r11     662     804   12574
r12      55     107    4551
r13      21      41     467
r14    1531    2928  292040
r15     141     384    5800
r16      55      61    1905
r17      61      92     848
r18      57      70    1782
r19     180     327   11801
r20     388     747   39758
r21      40      83     238
r22      64     114     783
r23     182     329    9550
r24     147     284    7310
r25     435     993    9091
r26     238     523   13346
r27     149     329    4017
r28      90     182    2420
r29      99     259    1546
r30     162     340    6901
r31     268     567   17210
r32     184     856  285600
r33     208     509   10887
r34     116     118    3387
r35     287     343    3894
r36     459     799   68143
r37     201     452    5146
r38      13      34      28
r39      21      53      51
r40     239     687    8340
r41      32     105     111
r42      94     185     805
r43      34      83     146
r44       7      20       5
r45     454     567   15975
r46     130     260    1388
r47     299     573   11248
r48     336     447   16576
r49       9      24      11
r50      84     117    1766
r51     355    1047  241508
r52     210     482   18004
r53      35      68     466
r54       8      14      22
r55      31      61     372
r56      79     168    2372
r57     843    1841  304803
r58      12      23     197
17

sum_mao[1,3]
         c1       c2       c3
r1    10358    16554  1115348

mao[58,3]
         c1      c2      c3
 r1     185     296   12466
 r2       3       6       1
 r3      27      50     466
 r4      55     100    2283
 r5     135     253   33388
 r6      19      34     212
 r7      19      36     215
 r8      10      20       8
 r9      64     126    7812
r10      21      34     197
r11     652     792   12187
r12      59     103    4607
r13      22      38     453
r14    1554    2573  174872
r15     150     213    5633
r16      56      72    1742
r17      61     100     748
r18      58      82    1627
r19     182     212   11456
r20     384     481   30163
r21      36      44      68
r22      65      75     573
r23     185     218    8371
r24     149     190    6324
r25     433    1010    9883
r26     237     391   10334
r27     144     270    3116
r28      93     131    2101
r29     102     221    1247
r30     167     262    5955
r31     267     380   12206
r32     184     856  285600
r33     207     544    9672
r34     119     127    3146
r35     289     346    3949
r36     456     682   38561
r37     204     260    2067
r38      14      26      17
r39      21      38      19
r40     235     539    5203
r41      35     100      92
r42      95     130     338
r43      33      53      32
r44       7      18       4
r45     452     638   69100
r46     130     253    1217
r47     296     548    9928
r48     336     447   16576
r49       8      15       4
r50      85     140    2892
r51     343     456   17176
r52     211     297   15581
r53      36      48     513
r54       9      12      33
r55      34      46     461
r56      83     103    2544
r57     829     996  269705
r58      13      23     204
18

sum_mao[1,3]
        c1      c2      c3
r1   10378   16430  840216

mao[58,3]
         c1      c2      c3
 r1     196     231   12490
 r2       2       4       3
 r3      27      46     586
 r4      62      93    2888
 r5     135     212   29611
 r6      18      31     235
 r7      18      32     294
 r8      10      19      83
 r9      59     111    2738
r10      21      37     141
r11     668     811   12652
r12      56      75    3696
r13      21      27     366
r14    1556    1747   41202
r15     151     252    4916
r16      55      69    1488
r17      60      97     640
r18      56      78    1343
r19     179     216    6534
r20     387     441    3412
r21      35      62     105
r22      64      85     314
r23     185     207    2625
r24     150     183    2887
r25     434    1015    9994
r26     234     404    6447
r27     146     282    2509
r28      88     122     908
r29     101     239    1126
r30     160     255    3386
r31     269     408    5284
r32     183     847  278880
r33     206     538    8296
r34     119     123    2788
r35     284     340    3812
r36     461     643   23515
r37     199     403    3363
r38      14      32      21
r39      21      49      37
r40     238     628    5858
r41      32     101      84
r42      95     173     583
r43      36      81     114
r44       7      20       4
r45     455     680   38363
r46     128     261    1054
r47     298     571   11168
r48     337     448   16675
r49      10      24      10
r50      86     146    1425
r51     349     578   44210
r52     203     401   26488
r53      35      59     232
r54       9      15      16
r55      33      56     201
r56      83     127    1639
r57     844    1182  210370
r58      10      13     107
19

sum_mao[1,3]
         c1       c2       c3
r1    10354    18795  1202281

mao[58,3]
         c1      c2      c3
 r1     194     307   15766
 r2       2       6       1
 r3      26      65     493
 r4      60     132    3121
 r5     134     257   37821
 r6      18      44     216
 r7      18      48     221
 r8       9      26       8
 r9      59     172   19608
r10      21      42     150
r11     656     797   12303
r12      58      98    5119
r13      17      29     311
r14    1530    2499  172913
r15     152     281    6392
r16      57      69    2080
r17      62      97     888
r18      56      74    1746
r19     177     273   11274
r20     388     693   40159
r21      37      75     207
r22      65     110     818
r23     186     305   10039
r24     152     261    7872
r25     430     982    9056
r26     238     498   13529
r27     149     318    4078
r28      90     167    2446
r29     101     261    1636
r30     166     327    7339
r31     271     543   17821
r32     178     826  265680
r33     201     507   10268
r34     116     122    3445
r35     288     345    3921
r36     463     806   69345
r37     204     451    5386
r38      14      36      33
r39      21      52      52
r40     241     685    8619
r41      34     111     127
r42      91     177     766
r43      36      87     166
r44       7      20       5
r45     451     563   15764
r46     129     260    1389
r47     297     569   11088
r48     336     447   16576
r49      10      26      14
r50      87     136    1814
r51     348     577   43955
r52     203     466   41037
r53      37      69     341
r54      12      21      38
r55      33      63     263
r56      83     143    2239
r57     843    1354  294319
r58      12      20     200
20

sum_mao[1,3]
         c1       c2       c3
r1    10359    18225  1176670

mao[58,3]
         c1      c2      c3
 r1     188     348   12126
 r2       3       6       8
 r3      27      56     623
 r4      60     123    2859
 r5     131     261   29409
 r6      19      39     278
 r7      18      38     315
 r8      10      21      92
 r9      60     124    3137
r10      21      38     201
r11     670     842    7655
r12      56     108    3899
r13      21      41     386
r14    1536    2937  293976
r15     143     256    4650
r16      56      91    1629
r17      62     126     726
r18      57     101    1473
r19     177     209    7079
r20     393     448    3547
r21      32      56      92
r22      64      84     334
r23     186     207    2867
r24     149     179    3008
r25     432     992    9264
r26     245     416    7469
r27     145     275    2613
r28      90     123    1001
r29     100     233    1167
r30     162     254    3665
r31     266     399    5462
r32     184     856  285600
r33     205     574    8978
r34     118     166    2891
r35     292     483   30826
r36     457     683   38732
r37     200     402    3589
r38      13      29      19
r39      21      49      39
r40     241     630    6355
r41      32     101      89
r42      96     174     628
r43      35      78     114
r44       6      17       3
r45     454     681   62211
r46     129     266    1138
r47     301     577   11409
r48     341     454   17075
r49      10      24      10
r50      86     150    2930
r51     345     572   43193
r52     206     407   27291
r53      34      57     263
r54      11      18      34
r55      31      53     214
r56      82     123    1746
r57     837    1145  220092
r58      13      25     191

sum_mao_mit[1,3]
           c1         c2         c3
r1      10368      18341  1828173.2
1352.0995

mao_mit[58,3]
            c1         c2         c3
 r1        190        316  19273.493
 r2          3          6  11.667895
 r3         28         57  901.86461
 r4         61        116  3900.7683
 r5        134        248  37734.137
 r6         19         38  368.75934
 r7         18         38  444.52618
 r8          9         20  104.61303
 r9         58        130  11178.496
r10         21         38  268.66263
r11        664        813  11683.991
r12         56         99  5227.2572
r13         20         33  437.81039
r14       1545       2494  524421.68
r15        151        303  11387.942
r16         56         78    2772.59
r17         62        113  6650.5447
r18         56         87  3006.6678
r19        179        263   13831.28
r20        388        596   53926.68
r21         36         59  420.21039
r22         65         95  1209.4711
r23        185        269  12383.815
r24        149        232  9084.3209
r25        433        996  9501.3099
r26        239        448  16869.929
r27        148        297  4605.3584
r28         90        150  2968.3205
r29        101        237  2068.4337
r30        161        288    8260.91
r31        269        463  24525.949
r32        182        847  279749.39
r33        204        549  14079.224
r34        118        135   4728.362
r35        287        363  10462.263
r36        459        702  45531.986
r37        201        353   15232.58
r38         14         30  54.776184
r39         21         45  113.83039
r40        239        611  12915.702
r41         33        101  183.51566
r42         94        160    1868.34
r43         34         69  416.84671
r44          7         20  7.6014474
r45        452        655  55716.535
r46        130        266  1707.5225
r47        297        564  11175.783
r48        339        479  25945.075
r49          9         22  34.876842
r50         86        145  2864.1316
r51        349        716  139370.94
r52        206        402  27597.749
r53         35         60  526.69158
r54         10         17  54.318289
r55         33         58  511.55197
r56         82        138  3064.4172
r57        840       1390  374557.23
r58         13         24  270.45408
11
5157
1

sum_mao[1,3]
        c1      c2      c3
r1     436    1212  120898

mao[58,3]
        c1     c2     c3
 r1     34     76    803
 r2      2      5     24
 r3      0      0      0
 r4      4     18     82
 r5      4      9     98
 r6      0      0      0
 r7      4     12     97
 r8      3      7     53
 r9      3      7     52
r10      0      0      0
r11     16     26    772
r12     49    116  14593
r13     15     49   1201
r14     62    407  70208
r15      5     14    140
r16      1      2      4
r17      0      0      0
r18      0      0      0
r19      0      0      0
r20      1      2      1
r21      1      1      1
r22      1      1      1
r23      1      2      2
r24      0      0      0
r25      3      5     16
r26      3      5      8
r27      0      0      0
r28      3      5      9
r29      0      0      0
r30      2      3      2
r31      1      2      1
r32     13     23     88
r33      7     12     96
r34      1      1      2
r35      5      7     30
r36     12     13      2
r37     11     14     29
r38      0      0      0
r39      1      1      0
r40      4      7     28
r41      3      5      4
r42      1      1      1
r43      4      5      2
r44      0      0      0
r45      4      7     28
r46      0      0      0
r47      3      7     55
r48      3      4      9
r49      2      3      1
r50      0      0      0
r51     69    163  26289
r52     21     49   2087
r53      2      4     14
r54      1      2      3
r55      1      2      4
r56     14     28    772
r57     20     39   1664
r58     16     41   1522
2

sum_mao[1,3]
        c1      c2      c3
r1     441    1220  132972

mao[58,3]
        c1     c2     c3
 r1     34     83   1040
 r2      2      5     28
 r3      1      2      7
 r4      3      8     62
 r5      4      9    111
 r6      0      0      0
 r7      4      9    111
 r8      4      9    111
 r9      3      7     62
r10      0      0      0
r11     21     42   2827
r12     50    118  17328
r13     16     38   1774
r14     62    407  70208
r15      5     12    173
r16      1      2      7
r17      1      2      5
r18      0      0      0
r19      1      2      4
r20      3      5      8
r21      2      3      4
r22      1      2      2
r23      1      2      2
r24      0      0      0
r25      3      5     43
r26      2      3     12
r27      0      0      0
r28      3      5     23
r29      0      0      0
r30      2      3     13
r31      1      1      1
r32     11     16     28
r33      8     15    271
r34      1      2      5
r35      5      7     65
r36     11     11      1
r37     10     15     53
r38      0      0      0
r39      1      1      0
r40      2      3      6
r41      3      4     11
r42      0      0      0
r43      5      7      9
r44      0      0      0
r45      5     10     60
r46      0      0      0
r47      3      7     62
r48      3      5     20
r49      2      3      3
r50      0      0      0
r51     66    156  30193
r52     20     47   2753
r53      2      5     25
r54      2      4     21
r55      2      5     25
r56     14     32   1271
r57     18     41   2121
r58     17     40   2003
3

sum_mao[1,3]
        c1      c2      c3
r1     451    2044  810222

mao[58,3]
         c1      c2      c3
 r1      32      74     782
 r2       2      10     103
 r3       0       0       0
 r4       3      30     190
 r5       3      10     290
 r6       0       0       0
 r7       4      29     432
 r8       4      20     408
 r9       5      25     619
r10       0       0       0
r11      18      30    3545
r12      51     208   74157
r13      16     101    5624
r14      63     814  612128
r15       7      33     976
r16       1       2      17
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       3       5       8
r21       1       1       2
r22       1       2       2
r23       1       2       2
r24       1       1       3
r25       3       5      56
r26       4       7      28
r27       0       0       0
r28       4       5       2
r29       0       0       0
r30       2       3       5
r31       1       2       2
r32      13      23      88
r33       7      12     342
r34       1       1       9
r35       5       7     109
r36      13      14       2
r37      12      15     117
r38       0       0       0
r39       1       1       1
r40       2       3       9
r41       3       4      16
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       7      16     160
r46       0       0       0
r47       4      11     502
r48       3       5      30
r49       2       3       5
r50       0       0       0
r51      66     245   85672
r52      20      80    6777
r53       2       6      49
r54       1       3      10
r55       2       6      50
r56      14      38    2575
r57      22      50    6607
r58      17      77    7709
4

sum_mao[1,3]
       c1     c2     c3
r1    431    995  34689

mao[58,3]
        c1     c2     c3
 r1     31     64    532
 r2      2      6      8
 r3      0      0      0
 r4      3     13     14
 r5      3      6     24
 r6      1      4      2
 r7      4     15     33
 r8      3      9     17
 r9      4     13     30
r10      0      0      0
r11     15     15    190
r12     51     69    251
r13     15     37    192
r14     52    257  26343
r15      5     14     34
r16      1      1      1
r17      0      0      0
r18      0      0      0
r19      0      0      0
r20      5      9     32
r21      1      1      0
r22      1      2      0
r23      1      2      0
r24      0      0      0
r25      2      2      2
r26      2      2      1
r27      1      1      0
r28      4      5      2
r29      0      0      0
r30      2      2      1
r31      1      1      0
r32     12     20     54
r33      8     13     32
r34      1      1      1
r35      6      8     11
r36     14     16      5
r37     11     14      7
r38      0      0      0
r39      1      1      0
r40      2      3      1
r41      3      4      1
r42      1      1      0
r43      4      5      2
r44      0      0      0
r45      5     10     60
r46      0      0      0
r47      2      3     11
r48      2      3      1
r49      2      2      0
r50      0      0      0
r51     64    161   5484
r52     21     56    512
r53      5     12     21
r54      1      2      1
r55      2      5      3
r56     13     29    153
r57     23     46    498
r58     18     30    122
5

sum_mao[1,3]
         c1       c2       c3
r1      448     1965  1058419

mao[58,3]
         c1      c2      c3
 r1      34      86    1142
 r2       2       6     154
 r3       0       0       0
 r4       2      18     134
 r5       3       9     362
 r6       0       0       0
 r7       5      25     993
 r8       3      10     344
 r9       5      16     942
r10       0       0       0
r11      14      23    4606
r12      49     165   95550
r13      15      84    7867
r14      66     898  748034
r15       6      25    1300
r16       1       2      29
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       2       3       2
r21       1       1       3
r22       1       1       4
r23       2       3       2
r24       1       1       7
r25       3       4     135
r26       2       3      34
r27       1       1      12
r28       3       4      68
r29       0       0       0
r30       3       4      90
r31       2       3      17
r32      13      23      88
r33       7      14     641
r34       1       2      16
r35       5       8     200
r36      12      13       2
r37      11      16     184
r38       1       1       4
r39       1       1       1
r40       2       3      19
r41       3       5      30
r42       0       0       0
r43       5       7       9
r44       0       0       0
r45       6      13     104
r46       0       0       0
r47       4      13     644
r48       2       3      25
r49       2       3       9
r50       0       0       0
r51      65     211  148411
r52      20      63   12078
r53       3       7     197
r54       1       2      18
r55       1       3      22
r56      13      30    3933
r57      23      50   12789
r58      21      79   17164
6

sum_mao[1,3]
        c1      c2      c3
r1     442    2101  935239

mao[58,3]
         c1      c2      c3
 r1      38     100    1578
 r2       2      14     119
 r3       0       0       0
 r4       2      25      96
 r5       4      16     645
 r6       0       0       0
 r7       5      48     779
 r8       3      21     265
 r9       5      34     709
r10       0       0       0
r11      13      15    2077
r12      50      59     107
r13      15      84    2921
r14      61     934  816316
r15       7      42    1013
r16       1       1      20
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       2       3       2
r21       1       1       2
r22       1       1       2
r23       1       1       3
r24       0       0       0
r25       3       4      61
r26       3       5       8
r27       1       1       4
r28       3       4      25
r29       0       0       0
r30       1       1       2
r31       2       3       7
r32      15      30     159
r33       7      12     367
r34       1       1      10
r35       9      12     380
r36      12      13       2
r37      11      14     105
r38       0       0       0
r39       1       1       1
r40       2       3      10
r41       3       5      17
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       5      10      60
r46       0       0       0
r47       3       8     371
r48       2       3      14
r49       2       3       5
r50       0       0       0
r51      66     308   89710
r52      21     103    7820
r53       4      15     206
r54       1       3      11
r55       1       4      13
r56      13      42    2369
r57      19      52    5383
r58      16      37    1463
7

sum_mao[1,3]
        c1      c2      c3
r1     450    1933  632716

mao[58,3]
         c1      c2      c3
 r1      32      70     683
 r2       2      13      55
 r3       0       0       0
 r4       2      22      44
 r5       6      27     718
 r6       0       0       0
 r7       5      43     359
 r8       3      19     122
 r9       3      19     117
r10       0       0       0
r11      15      40    1040
r12      54      90     648
r13      15      78    1298
r14      60     790  577490
r15       6      31     326
r16       1       1       9
r17       1       1       4
r18       0       0       0
r19       0       0       0
r20       2       3       2
r21       1       1       1
r22       1       1       1
r23       1       1       1
r24       0       0       0
r25       2       4       7
r26       3       5       8
r27       1       2       1
r28       4       5       2
r29       0       0       0
r30       2       3       3
r31       1       1       1
r32      13      23      88
r33       7      12     162
r34       1       1       4
r35       6       8      75
r36      11      11       1
r37      11      14      47
r38       0       0       0
r39       1       1       0
r40       4       7      28
r41       3       5       7
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       6      14     126
r46       0       0       0
r47       6      14     686
r48       3       4      14
r49       2       3       2
r50       0       0       0
r51      67     280   40194
r52      20      90    3120
r53       3      10      51
r54       1       3       5
r55       2       7      23
r56      15      45    1360
r57      21      55    2764
r58      20      51    1017
8

sum_mao[1,3]
        c1      c2      c3
r1     458    1946  670635

mao[58,3]
         c1      c2      c3
 r1      35      87    1129
 r2       2      14      54
 r3       1      10      13
 r4       4      46     172
 r5       5      24     511
 r6       0       0       0
 r7       4      36     226
 r8       4      27     213
 r9       3      20     115
r10       0       0       0
r11      16      18    1341
r12      51      60     112
r13      16      81    1418
r14      60     815  616140
r15       9      48     698
r16       2       3       2
r17       1       1       3
r18       1       1       2
r19       0       0       0
r20       1       1       1
r21       1       1       1
r22       1       1       1
r23       2       3       2
r24       1       1       2
r25       3       4      25
r26       3       5       8
r27       0       0       0
r28       3       4      11
r29       0       0       0
r30       2       3       3
r31       1       1       1
r32      11      16      28
r33       7      11     147
r34       1       1       3
r35       6       7      60
r36      13      14       2
r37      11      14      45
r38       0       0       0
r39       1       1       0
r40       3       5       8
r41       3       4       6
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       6      13     104
r46       0       0       0
r47       4       3     318
r48       5       7      38
r49       2       3       2
r50       0       0       0
r51      68     282   39455
r52      20      90    2956
r53       2       7      21
r54       1       3       4
r55       1       4       5
r56      14      41    1103
r57      23      56    3147
r58      20      44     977
9

sum_mao[1,3]
        c1      c2      c3
r1     454    1385  135231

mao[58,3]
         c1      c2      c3
 r1      34      80     926
 r2       2       8      20
 r3       0       0       0
 r4       3      19      36
 r5       7      19     328
 r6       0       0       0
 r7       4      20      83
 r8       4      16      78
 r9       5      20     118
r10       0       0       0
r11      17      22     691
r12      52      79     430
r13      15      52     485
r14      67     513  114639
r15       7      25     169
r16       1       1       3
r17       0       0       0
r18       1       1       3
r19       0       0       0
r20       2       3       2
r21       1       1       0
r22       1       1       0
r23       1       1       0
r24       0       0       0
r25       3       4      12
r26       2       3       3
r27       0       0       0
r28       3       4       7
r29       0       0       0
r30       1       1       1
r31       5       7       9
r32      12      20      54
r33       7      10      61
r34       1       1       2
r35       5       6      19
r36      13      15       3
r37      11      14      17
r38       0       0       0
r39       1       1       0
r40       4       7      28
r41       3       5       2
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       4       7      28
r46       1       1       1
r47       5      10     172
r48       2       3       2
r49       2       3       1
r50       0       0       0
r51      63     190   13645
r52      20      62    1179
r53       2       5       9
r54       1       2       2
r55       1       3       2
r56      13      30     396
r57      22      46    1220
r58      19      39     343
10

sum_mao[1,3]
        c1      c2      c3
r1     444    1841  562577

mao[58,3]
         c1      c2      c3
 r1      35      87    1129
 r2       2      10      84
 r3       0       0       0
 r4       2      19      68
 r5       4      12     417
 r6       0       0       0
 r7       4      27     350
 r8       3      14     186
 r9       4      19     321
r10       0       0       0
r11      18      29    2871
r12      51     197   60042
r13      16      96    4552
r14      54     661  401888
r15       6      27     575
r16       1       2      14
r17       0       0       0
r18       1       2      13
r19       0       0       0
r20       2       3       2
r21       1       1       1
r22       1       1       2
r23       2       3       2
r24       1       1       2
r25       3       5      46
r26       4       7      28
r27       0       0       0
r28       4       5       2
r29       0       0       0
r30       3       5      10
r31       2       3       5
r32      13      22      66
r33       8      13     361
r34       1       1       7
r35       5       7      88
r36      10      10       0
r37      13      16     111
r38       0       0       0
r39       1       1       1
r40       2       3       7
r41       3       4      13
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       7      16     160
r46       0       0       0
r47       4      10     407
r48       3       4      24
r49       2       3       4
r50       0       0       0
r51      67     239   70454
r52      21      80    5968
r53       3       9      87
r54       3       9      71
r55       1       3      10
r56      13      33    1720
r57      18      39    3409
r58      18      78    6997
11

sum_mao[1,3]
        c1      c2      c3
r1     452    1943  769591

mao[58,3]
         c1      c2      c3
 r1      33      80     954
 r2       2      11     128
 r3       2      18     123
 r4       2      22     106
 r5       6      18    1399
 r6       0       0       0
 r7       5      38     835
 r8       3      15     285
 r9       3      15     277
r10       0       0       0
r11      14      37     888
r12      51      60     112
r13      15      72    3471
r14      60     815  616140
r15       7      38    1226
r16       1       2      22
r17       2       3      44
r18       1       2      20
r19       0       0       0
r20       1       2       8
r21       3       4      24
r22       1       1       5
r23       1       2       9
r24       0       0       0
r25       2       4      22
r26       2       3      29
r27       0       0       0
r28       3       5      83
r29       0       0       0
r30       2       4      33
r31       4       6     104
r32      12      20      54
r33       7      10     432
r34       1       2      11
r35       6       8     197
r36      12      13       2
r37      12      15     155
r38       0       0       0
r39       1       1       1
r40       2       3      12
r41       3       4      20
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       4       7      28
r46       1       1       4
r47       6      20    1373
r48       3       4      38
r49       2       3       6
r50       0       0       0
r51      66     286  110756
r52      26     111   14686
r53       2       6      64
r54       2       6      51
r55       1       3      16
r56      14      42    3627
r57      22      60    9761
r58      17      36    1948
12

sum_mao[1,3]
        c1      c2      c3
r1     467    1338  157608

mao[58,3]
         c1      c2      c3
 r1      36      93    1351
 r2       2       5      36
 r3       0       0       0
 r4       4      21     125
 r5       4      10     151
 r6       0       0       0
 r7       6      20     335
 r8       3       8      81
 r9       6      15     318
r10       0       0       0
r11      16      24    1546
r12      54      73     283
r13      16      49    1416
r14      65     492  105272
r15      11      35    1038
r16       1       2       7
r17       0       0       0
r18       0       0       0
r19       1       2       5
r20       1       1       3
r21       1       1       1
r22       1       1       2
r23       1       2       3
r24       1       2       4
r25       2       3      16
r26       2       3      13
r27       0       0       0
r28       3       4      36
r29       0       0       0
r30       2       3      16
r31       1       1       3
r32      12      20      54
r33       7      12     153
r34       1       1       4
r35       7      10      92
r36      12      13       2
r37      15      19      85
r38       0       0       0
r39       1       1       0
r40       2       3       5
r41       3       4       7
r42       2       3       3
r43       4       5       2
r44       0       0       0
r45       5      10      60
r46       1       2       2
r47       3       8      85
r48       7      11      73
r49       2       2       2
r50       0       0       0
r51      66     177   37300
r52      19      50    2635
r53       3       7      49
r54       1       2       4
r55       1       2       6
r56      14      30    1249
r57      20      42    2808
r58      19      34     867
13

sum_mao[1,3]
        c1      c2      c3
r1     448    2072  880105

mao[58,3]
         c1      c2      c3
 r1      36      90    1222
 r2       2      15      84
 r3       0       0       0
 r4       3      37     152
 r5       5      23     764
 r6       0       0       0
 r7       6      57     794
 r8       3      21     187
 r9       3      21     180
r10       0       0       0
r11      17      14    2508
r12      52      70     261
r13      16      89    2261
r14      62     922  793842
r15       9      51    1122
r16       2       3       2
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       2       3       2
r21       1       1       1
r22       1       1       2
r23       1       1       2
r24       0       0       0
r25       2       2      21
r26       2       2       8
r27       0       0       0
r28       4       5       2
r29       0       0       0
r30       3       3      16
r31       1       1       1
r32      11      16      28
r33       7      11     235
r34       1       1       4
r35       6       7      96
r36      11      11       1
r37      11      14      72
r38       0       0       0
r39       1       1       0
r40       2       2       7
r41       3       4      12
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       5      10      60
r46       0       0       0
r47       3       5     268
r48       4       6      39
r49       2       2       4
r50       0       0       0
r51      68     298   63449
r52      21      99    5258
r53       2       7      35
r54       1       3       7
r55       1       4       9
r56      14      42    1816
r57      21      53    4273
r58      16      39     996
14

sum_mao[1,3]
        c1      c2      c3
r1     432    1387  213537

mao[58,3]
         c1      c2      c3
 r1      32      74     782
 r2       2      10      30
 r3       0       0       0
 r4       4      32      95
 r5       5      16     269
 r6       0       0       0
 r7       5      32     194
 r8       4      20     117
 r9       5      25     176
r10       0       0       0
r11      16      20     744
r12      51      60     114
r13      15      57     700
r14      48     454  184778
r15       5      21     122
r16       1       1       5
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       3       5       8
r21       1       1       0
r22       1       2       1
r23       1       2       1
r24       0       0       0
r25       3       4      14
r26       4       7      28
r27       0       0       0
r28       4       5       2
r29       0       0       0
r30       2       3       1
r31       1       2       0
r32      14      27     130
r33       7      13      88
r34       1       1       2
r35       7       9      55
r36      13      15       3
r37      12      15      30
r38       0       0       0
r39       1       1       0
r40       3       5       8
r41       3       5       4
r42       1       1       0
r43       4       5       2
r44       0       0       0
r45       6      13     104
r46       0       0       0
r47       4       5     178
r48       2       3       3
r49       2       3       1
r50       1       2       2
r51      66     223   21106
r52      19      69    1518
r53       2       6      12
r54       1       3       3
r55       1       3       3
r56      13      33     552
r57      18      40    1108
r58      18      34     444
15

sum_mao[1,3]
        c1      c2      c3
r1     439    1577  306305

mao[58,3]
         c1      c2      c3
 r1      39     104    1698
 r2       2      12      39
 r3       0       0       0
 r4       3      27      70
 r5       3      10     127
 r6       0       0       0
 r7       4      30     163
 r8       3      18      87
 r9       5      29     231
r10       0       0       0
r11      14      14     749
r12      50      59     107
r13      15      63     918
r14      51     544  268736
r15       9      43     516
r16       1       1       6
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       3       5       8
r21       1       1       1
r22       1       2       1
r23       2       3       2
r24       2       3       3
r25       3       3      19
r26       3       5       8
r27       0       0       0
r28       4       5       2
r29       0       0       0
r30       2       3       2
r31       1       1       1
r32      12      18      35
r33       7      12     116
r34       2       2      12
r35       6       7      53
r36      11      11       1
r37      10      13      27
r38       0       0       0
r39       1       1       0
r40       2       3       3
r41       3       4       5
r42       1       1       1
r43       4       5       2
r44       0       0       0
r45       9      26     468
r46       0       0       0
r47       5       6     364
r48       2       3       4
r49       2       2       2
r50       0       0       0
r51      65     251   26514
r52      20      86    2180
r53       2       7      16
r54       1       3       3
r55       2       7      16
r56      15      44     925
r57      19      46    1545
r58      17      34     519
16

sum_mao[1,3]
        c1      c2      c3
r1     454    1722  420713

mao[58,3]
         c1      c2      c3
 r1      35      87    1129
 r2       2      12      60
 r3       0       0       0
 r4       3      28     110
 r5       4      13     327
 r6       0       0       0
 r7       5      38     395
 r8       3      18     135
 r9       4      24     230
r10       0       0       0
r11      17      19    1801
r12      51      60     112
r13      15      64    1483
r14      57     628  359216
r15      10      50    1039
r16       1       1      10
r17       2       3      19
r18       0       0       0
r19       0       0       0
r20       1       1       1
r21       1       1       1
r22       1       2       1
r23       2       3       2
r24       0       0       0
r25       2       2      14
r26       3       5       8
r27       1       1       2
r28       3       4      13
r29       1       1       2
r30       3       4       9
r31       1       1       1
r32      12      20      54
r33       7      13     186
r34       2       2      20
r35       8      10     152
r36      10      10       0
r37      12      18      64
r38       0       0       0
r39       1       1       0
r40       2       3       5
r41       3       5       9
r42       0       0       0
r43       5       7       9
r44       0       0       0
r45      10      30     630
r46       0       0       0
r47       5      11     523
r48       2       4       7
r49       2       3       3
r50       0       0       0
r51      65     263   43403
r52      23     102    4691
r53       2       7      25
r54       2       7      21
r55       2       7      26
r56      15      46    1519
r57      18      47    2306
r58      18      36     940
17

sum_mao[1,3]
        c1      c2      c3
r1     456    1647  349273

mao[58,3]
         c1      c2      c3
 r1      36      93    1351
 r2       2      11      52
 r3       0       0       0
 r4       4      34     169
 r5       5      16     442
 r6       0       0       0
 r7       5      34     342
 r8       3      15     116
 r9       3      15     112
r10       0       0       0
r11      18      27    1746
r12      54      73     283
r13      15      63    1281
r14      55     570  294120
r15       5      24     229
r16       1       1       9
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       1       2       2
r21       3       5       8
r22       1       2       2
r23       1       2       3
r24       3       6      24
r25       2       3      12
r26       4       7      28
r27       0       0       0
r28       3       6      14
r29       0       0       0
r30       1       2       1
r31       3       5      11
r32      12      18      35
r33       7      11     161
r34       2       2      17
r35       9      11     166
r36      13      15       3
r37      11      16      47
r38       0       0       0
r39       1       1       0
r40       2       2       3
r41       3       4       7
r42       0       0       0
r43       5       7       9
r44       0       0       0
r45       6      13     104
r46       0       0       0
r47       5      11     452
r48       3       5      14
r49       3       3       1
r50       0       0       0
r51      65     260   38822
r52      21      86    3477
r53       3      10      52
r54       2       6      19
r55       2       7      24
r56      13      41    1093
r57      23      66    3686
r58      17      36     724
18

sum_mao[1,3]
        c1      c2      c3
r1     441    1507  242537

mao[58,3]
         c1      c2      c3
 r1      34      86    1142
 r2       2      12      30
 r3       0       0       0
 r4       2      18      24
 r5       4      16     181
 r6       0       0       0
 r7       4      29     125
 r8       4      24     118
 r9       4      23     113
r10       0       0       0
r11      14      16     568
r12      53      71     263
r13      17      71     886
r14      48     488  215788
r15       6      27     172
r16       2       3       2
r17       0       0       0
r18       0       0       0
r19       1       2       1
r20       1       1       0
r21       2       3       2
r22       1       1       1
r23       2       3       2
r24       1       1       1
r25       2       3       6
r26       3       5       8
r27       1       1       1
r28       3       4       6
r29       0       0       0
r30       1       1       0
r31       1       1       1
r32      13      23      88
r33       7      14      81
r34       1       1       1
r35       9      12      75
r36      14      16       5
r37      15      19      47
r38       0       0       0
r39       1       1       0
r40       4       7      28
r41       4       7       6
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       9      26     468
r46       0       0       0
r47       5       4     275
r48       4       7      13
r49       2       3       1
r50       0       0       0
r51      62     231   18153
r52      21      87    1804
r53       2       7      12
r54       2       7      10
r55       1       4       3
r56      14      40     611
r57      18      43    1067
r58      16      33     346
19

sum_mao[1,3]
        c1      c2      c3
r1     444    1440  196423

mao[58,3]
         c1      c2      c3
 r1      34      86    1142
 r2       2      10      35
 r3       1       6       8
 r4       3      22      64
 r5       9      25     964
 r6       0       0       0
 r7       4      24     147
 r8       3      15      78
 r9       5      25     210
r10       0       0       0
r11      15      17     817
r12      51      60     112
r13      15      52     863
r14      50     424  159520
r15       9      39     495
r16       1       1       6
r17       1       2       3
r18       0       0       0
r19       0       0       0
r20       3       5       8
r21       2       3       2
r22       1       2       1
r23       1       2       1
r24       0       0       0
r25       2       2       8
r26       3       5       8
r27       1       1       1
r28       3       5       7
r29       0       0       0
r30       2       3       2
r31       2       3       2
r32      14      27     130
r33       9      18     179
r34       1       1       3
r35       5       7      35
r36      10      10       0
r37      10      13      26
r38       0       0       0
r39       1       1       0
r40       2       3       3
r41       4       6       9
r42       1       1       1
r43       4       5       2
r44       1       2       1
r45       9      26     468
r46       1       2       1
r47       3       6     110
r48       2       3       4
r49       2       3       2
r50       0       0       0
r51      65     241   25702
r52      21      84    2310
r53       2       7      15
r54       1       3       3
r55       2       7      16
r56      13      40     700
r57      19      51    1590
r58      19      34     609
20

sum_mao[1,3]
        c1      c2      c3
r1     427    1455  267511

mao[58,3]
         c1      c2      c3
 r1      31      68     627
 r2       2      12      38
 r3       1       8       9
 r4       2      17      31
 r5       5      17     346
 r6       0       0       0
 r7       4      28     160
 r8       3      17      85
 r9       3      17      81
r10       0       0       0
r11      15      12    1027
r12      51      60     112
r13      15      60     899
r14      50     505  230280
r15      10      44     632
r16       2       3       2
r17       0       0       0
r18       0       0       0
r19       0       0       0
r20       2       3       2
r21       1       2       0
r22       1       2       1
r23       1       2       1
r24       0       0       0
r25       2       2      11
r26       2       2       6
r27       1       1       2
r28       3       4      13
r29       0       0       0
r30       3       3      16
r31       1       1       1
r32      11      16      28
r33       8      15     140
r34       1       1       2
r35       5       6      30
r36      12      13       2
r37      11      28     504
r38       0       0       0
r39       1       1       0
r40       2       2       3
r41       4       5       9
r42       0       0       0
r43       4       5       2
r44       0       0       0
r45       7      18     216
r46       0       0       0
r47       5       6     339
r48       2       4       4
r49       2       3       2
r50       0       0       0
r51      65     231   26701
r52      20      79    2171
r53       2       7      16
r54       1       3       3
r55       2       7      16
r56      13      36     729
r57      18      41    1508
r58      20      38     704

sum_mao_mit[1,3]
           c1         c2         c3
r1        443       1636  493759.34
702.68011

mao_mit[58,3]
            c1         c2         c3
 r1         34         83  1169.5495
 r2          2         10  69.823553
 r3          0          2  31.829474
 r4          3         24  174.16632
 r5          5         15  462.25987
 r6          0          0        .94
 r7          5         30  495.21368
 r8          3         16  185.83039
 r9          4         19  298.22882
r10          0          0          0
r11         16         23  1700.8158
r12         51         90  15552.138
r13         15         67  2431.7263
r14         58        617  406411.29
r15          7         32   753.5225
r16          1          2  9.4993421
r17          0          1  5.0494737
r18          0          0  2.3531579
r19          0          0  1.0636842
r20          2          3  9.1756579
r21          1          2  4.3084211
r22          1          1  1.8735526
r23          1          2  2.6415789
r24          1          1  4.5768421
r25          3          4  28.571053
r26          3          4  17.247237
r27          0          0  1.5340789
r28          3          5  16.811974
r29          0          0      .1525
r30          2          3  12.435658
r31          2          2  11.887368
r32         12         21  85.260395
r33          7         13  216.33829
r34          1          1  6.9821053
r35          6          8  103.00316
r36         12         13  6.0698684
r37         12         16  100.49263
r38          0          0      .2525
r39          1          1        .25
r40          3          4  14.396184
r41          3          5  10.443553
r42          0          0  .89684211
r43          4          5  4.1073684
r44          0          0        .26
r45          6         15  226.29145
r46          0          0  .85315789
r47          4          8  378.25211
r48          3          5  22.941974
r49          2          3  2.9768421
r50          0          0        .31
r51         66        235  50335.953
r52         21         79  4680.1251
r53          3          8  55.705132
r54          1          4  17.876842
r55          1          5  18.167237
r56         14         37  1460.6653
r57         20         48   3519.983
r58         18         44  2623.9711
10
181

. 
. 
. *************************************************************************************
. *       Final reports
. *************************************************************************************
. 
. matrix mao_mit_tot=mao_mit_est+mao_mit_sl+mao_mit_ot

. matrix list mao_mit_tot

mao_mit_tot[58,3]
            c1         c2         c3
 r1        657       1500  41583.684
 r2         48         65  106.54658
 r3         66        107  969.33355
 r4        143        332  8050.1551
 r5        235        558  46252.883
 r6         43         66  378.22776
 r7         72        201  10300.586
 r8         39         66  298.64658
 r9        274       1385  891605.86
r10         40         58  270.66263
r11       1480       2942  80842.748
r12        296       1027  107338.07
r13         57        172  4899.7925
r14       3220       6875  999350.32
r15        264        664  20679.088
r16         87        231  20483.714
r17        117        230  11249.349
r18         77        120  3132.3993
r19        363       1047   237016.4
r20        777       1640  100781.16
r21        112        174  1684.8788
r22        136        169  1218.0451
r23        399        672  15542.467
r24        212        304  9214.9762
r25       1410       2468  13272.287
r26        426        795  20379.639
r27        240        420  5090.2946
r28        254       1783  2291265.7
r29        151        319  2309.8446
r30        264        426  8678.9357
r31        448       1049  52246.144
r32        597       1558  286810.65
r33        312        780  26002.933
r34        184        277  6073.9593
r35        785       1205  19349.524
r36       1124       2274  115650.62
r37        733       1238  23930.289
r38         31         64  251.09447
r39         55         80  116.08039
r40        392        810  13116.508
r41        145        412  25172.278
r42        156        234   1915.813
r43        121        244  4100.2376
r44         51        115   1208.183
r45        891       1527  69525.723
r46        204        485  6284.7251
r47        430        800  12546.779
r48        516        749   26572.94
r49         36         88  854.85882
r50        123        207  4111.9237
r51        572       1157   190163.8
r52        448       1400  232124.32
r53        143        194  721.70776
r54         33         51  108.52553
r55         88        128  585.85013
r56        171        288  5107.7682
r57       1081       3074  842945.93
r58        125        545  172874.81

. matrix sum_mit_tot=sum_mit_est+sum_mit_sl+sum_mit_ot

. *matrix list sum_mit_tot
. matrix sum_mit_all=sum_mit_est\sum_mit_sl\sum_mit_ot\sum_mit_tot

. matrix list sum_mit_all

sum_mit_all[4,3]
           c1         c2         c3
r1      11143      27872  4772788.2
r1      10368      18341  1828173.2
r1        443       1636  493759.34
r1      21954      47849  7094720.7

. matrix modelmat_m=modelmat_m_e,modelmat_m_s,modelmat_m_o

. matrix list modelmat_m

modelmat_m[7,3]
       c1     c1     c1
r1   3.85   6.45   9.25
r2      2      0      0
r3     .1    .05      0
r4  20.75   2.35     .1
r5   20.1      0      0
r6      1      0      0
r7   10.2   2.05      0

. matrix ones=J(1,7,1)

. matrix modelmatt=ones*modelmat_m

. matrix list modelmatt

modelmatt[1,3]
      c1    c1    c1
r1    58  10.9  9.35

. di diestim_mi_e " " diestim_mi_s " " diestim_mi_o
58 10.9 9.35

. di nestim_mi_e " " nestim_mi_s " " nestim_mi_o
11139.55 5064.6 174.9

. 
. *************************************************************************************
. *       Totals
. *************************************************************************************
. 
. /* Report Total */
. matrix mao=sum_mit_all

. 
. matrix smao=mao[1...,3]

. matrix bmao=mao[1...,2]

. matrix nmao=mao[1...,1]

. 
. *matrix tmao=bmao
. 
. matrix cil_mao=J(4,1,0)

. matrix ciu_mao=J(4,1,0)

. forval i= 1/4{
  2. matrix smao[`i',1]=sqrt(smao[`i',1])
  3. *matrix tmao[`i',1]=bmao[`i',1]/smao[`i',1]
. matrix cil_mao[`i',1]=int(bmao[`i',1]-signi*smao[`i',1]+0.5)
  4. matrix ciu_mao[`i',1]=int(bmao[`i',1]+signi*smao[`i',1]+0.5)
  5. }

. 
. 
. matrix mao=nmao,cil_mao,bmao,ciu_mao,bmao,smao

. matrix list mao

mao[4,6]
           c1         c1         c2         c1         c2         c3
r1      11143      23503      27872      32241      27872  2184.6712
r1      10368      15637      18341      21045      18341  1352.0995
r1        443        231       1636       3041       1636  702.68011
r1      21954      42522      47849      53176      47849  2663.5917

. 
. svmat mao, names(v)
number of observations will be reset to 4
Press any key to continue, or Break to abort
obs was 0, now 4

. ren v1 n

. ren v2 cl

. ren v3 e

. ren v4 cu

. ren v5 em

. ren v6 s

. 
. export excel n cl e cu em s using "$excel/${name}smi.xls", sheet("totals") sheetmodify firstrow(varlabels) missing(".")
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. clear

. 
. 
. /* Report Strata */
. foreach k in mao_mit_est mao_mit_sl mao_mit_ot mao_mit_tot{
  2. matrix mao=`k'
  3. 
. matrix smao=mao[1...,3]
  4. matrix bmao=mao[1...,2]
  5. matrix nmao=mao[1...,1]
  6. 
. *matrix tmao=bmao
. 
. matrix cil_mao=J(nstrata,1,0)
  7. matrix ciu_mao=J(nstrata,1,0)
  8. forval i= 1/58{
  9. matrix smao[`i',1]=sqrt(smao[`i',1])
 10. *matrix tmao[`i',1]=bmao[`i',1]/smao[`i',1]
. matrix cil_mao[`i',1]=int(bmao[`i',1]-signi*smao[`i',1]+0.5)
 11. matrix ciu_mao[`i',1]=int(bmao[`i',1]+signi*smao[`i',1]+0.5)
 12. }
 13. 
. 
. matrix mao=nmao,cil_mao,bmao,ciu_mao,bmao,smao
 14. matrix list mao
 15. 
. svmat mao, names(v)
 16. ren v1 n
 17. ren v2 cl
 18. ren v3 e
 19. ren v4 cu
 20. ren v5 em
 21. ren v6 s
 22. 
. export excel n cl e  cu em s using "$excel/${name}smi.xls", sheet("s`k'") sheetmodify firstrow(varlabels) missing(".")
 23. clear
 24. 
. 
. }

mao[58,6]
            c1         c1         c2         c1         c2         c3
 r1        433        810       1101       1392       1101  145.39822
 r2         43         39         49         59         49  5.0055101
 r3         38         36         48         60         48  5.9698805
 r4         79         66        192        318        192   63.04935
 r5         96        115        295        475        295  89.757928
 r6         24         22         28         34         28  2.9203461
 r7         49        -60        133        327        133  96.751466
 r8         27         24         30         36         30  2.8641156
 r9        212       -639       1236       3112       1236  938.15198
r10         19         17         20         23         20  1.4142136
r11        800       1587       2106       2625       2106  259.72667
r12        189        250        838       1426        838  294.20855
r13         22        -17         72        162         72   45.05836
r14       1617       3240       3764       4288       3764  261.75819
r15        106        144        329        514        329   92.39926
r16         30       -114        151        417        151  133.04745
r17         55        -19        116        252        116  67.777244
r18         21         11         33         55         33  11.107584
r19        184       -160        784       1729        784   472.4236
r20        387        608       1041       1474       1041  216.43776
r21         75         42        113        184        113  35.501549
r22         70         68         73         78         73  2.5885375
r23        213        289        401        513        401  56.178377
r24         62         49         71         93         71  11.228465
r25        974       1346       1468       1590       1468  61.175206
r26        184        225        343        461        343  59.097064
r27         92         79        123        167        123  21.986407
r28        161      -1396       1628       4653       1628  1512.7064
r29         50         51         82        113         82  15.532496
r30        101         95        135        175        135  20.139265
r31        177        251        584        917        584  166.45813
r32        403        523        690        857        690  83.522405
r33        101          2        218        434        218   108.2006
r34         65         68        141        214        141  36.587091
r35        492        647        834       1021        834  93.724371
r36        653       1029       1559       2089       1559  264.78776
r37        520        684        869       1054        869  92.721175
r38         17          6         34         62         34  14.002349
r39         33         31         34         37         34  1.4142136
r40        150        168        195        222        195  13.653181
r41        109         -9        306        622        306  158.04531
r42         62         60         74         88         74  6.8246747
r43         83         49        170        291        170  60.657098
r44         44         26         95        164         95  34.645657
r45        433        624        857       1090        857  116.54569
r46         74         84        219        354        219  67.648721
r47        129        165        228        291        228  31.507833
r48        174        216        265        314        265  24.595183
r49         25          6         63        120         63  28.583302
r50         37         -8         62        133         62  35.319713
r51        157        163        206        249        206  21.375232
r52        221         25        919       1813        919  447.04188
r53        105        102        126        150        126   11.80301
r54         22         18         30         42         30    6.02747
r55         54         50         65         80         65  7.4920572
r56         75         65        113        161        113  24.138883
r57        221        272       1636       3000       1636  681.81282
r58         94       -347        477       1302        477  412.28678
number of observations will be reset to 58
Press any key to continue, or Break to abort
obs was 0, now 58
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[58,6]
            c1         c1         c2         c1         c2         c3
 r1        190         38        316        594        316  138.82901
 r2          3          0          6         13          6    3.41583
 r3         28         -2         57        117         57  30.031061
 r4         61         -8        116        241        116  62.456131
 r5        134       -140        248        637        248  194.25277
 r6         19          0         38         76         38  19.203108
 r7         18         -3         38         80         38   21.08379
 r8          9          0         20         40         20  10.228051
 r9         58        -80        130        341        130  105.72841
r10         21          5         38         71         38  16.390931
r11        664        597        813       1029        813  108.09251
r12         56        -45         99        244         99  72.299773
r13         20         -8         33         75         33  20.923919
r14       1545       1046       2494       3942       2494  724.16965
r15        151         90        303        516        303   106.7143
r16         56        -26         78        183         78  52.655389
r17         62        -49        113        276        113  81.550872
r18         56        -22         87        197         87   54.83309
r19        179         28        263        498        263  117.60646
r20        388        132        596       1060        596  232.22119
r21         36         18         59        100         59  20.499034
r22         65         25         95        165         95   34.77745
r23        185         46        269        492        269  111.28259
r24        149         41        232        423        232  95.311704
r25        433        801        996       1191        996  97.474663
r26        239        188        448        708        448  129.88429
r27        148        161        297        433        297  67.862791
r28         90         41        150        259        150  54.482296
r29        101        146        237        328        237  45.480036
r30        161        106        288        470        288  90.889548
r31        269        150        463        776        463  156.60763
r32        182       -210        847       1905        847  528.91341
r33        204        312        549        786        549   118.6559
r34        118         -2        135        273        135  68.763086
r35        287        158        363        568        363   102.2852
r36        459        275        702       1129        702  213.38225
r37        201        106        353        600        353  123.42034
r38         14         15         30         45         30  7.4010934
r39         21         24         45         66         45  10.669133
r40        239        384        611        838        611  113.64727
r41         33         74        101        128        101  13.546795
r42         94         74        160        246        160  43.224299
r43         34         28         69        110         69  20.416824
r44          7         14         20         26         20  2.7570722
r45        452        183        655       1127        655   236.0435
r46        130        183        266        349        266  41.322179
r47        297        353        564        775        564  105.71558
r48        339        157        479        801        479  161.07475
r49          9         10         22         34         22  5.9056619
r50         86         38        145        252        145  53.517582
r51        349        -30        716       1463        716  373.32418
r52        206         70        402        734        402   166.1257
r53         35         14         60        106         60  22.949762
r54         10          2         17         32         17  7.3700943
r55         33         13         58        103         58  22.617515
r56         82         27        138        249        138  55.357179
r57        840        166       1390       2614       1390  612.01081
r58         13         -8         24         57         24  16.445488
number of observations will be reset to 58
Press any key to continue, or Break to abort
obs was 0, now 58
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[58,6]
            c1         c1         c2         c1         c2         c3
 r1         34         15         83        151         83  34.198676
 r2          2         -6         10         27         10  8.3560489
 r3          0         -8          2         13          2  5.6417616
 r4          3         -1         24         50         24  13.197209
 r5          5        -27         15         58         15  21.500229
 r6          0         -1          0          2          0  .96953597
 r7          5        -14         30         75         30  22.253397
 r8          3        -10         16         43         16  13.631962
 r9          4        -15         19         54         19  17.269303
r10          0          0          0          0          0          0
r11         16        -58         23        105         23  41.240948
r12         51       -158         90        339         90  124.70821
r13         15        -31         67        166         67  49.312537
r14         58       -657        617       1892        617  637.50395
r15          7        -22         32         87         32  27.450364
r16          1         -3          2          8          2  3.0821003
r17          0         -2          1          5          1  2.2471034
r18          0         -2          0          3          0  1.5340006
r19          0         -1          0          2          0  1.0313507
r20          2         -2          3          9          3  3.0291348
r21          1         -1          2          6          2  2.0756736
r22          1         -1          1          4          1  1.3687778
r23          1          0          2          5          2  1.6252935
r24          1         -2          1          5          1  2.1393555
r25          3         -6          4         15          4  5.3451897
r26          3         -3          4         12          4  4.1529793
r27          0         -1          0          2          0  1.2385794
r28          3         -2          5         13          5  4.1002407
r29          0          0          0          1          0  .39051248
r30          2         -3          3         10          3  3.5264228
r31          2         -4          2          9          2  3.4478063
r32         12          3         21         39         21  9.2336555
r33          7        -15         13         42         13  14.708443
r34          1         -3          1          6          1  2.6423674
r35          6        -11          8         28          8  10.149047
r36         12          8         13         18         13  2.4637103
r37         12         -3         16         36         16  10.024601
r38          0          0          0          1          0  .50249378
r39          1          0          1          2          1         .5
r40          3         -3          4         12          4  3.7942304
r41          3          0          5         11          5  3.2316486
r42          0         -1          0          2          0  .94701748
r43          4          1          5          9          5  2.0266644
r44          0          0          0          1          0  .50990195
r45          6        -14         15         45         15  15.042987
r46          0         -1          0          2          0  .92366547
r47          4        -30          8         47          8  19.448704
r48          3         -4          5         15          5   4.789778
r49          2          0          3          6          3  1.7253527
r50          0          0          0          1          0  .55677644
r51         66       -213        235        684        235  224.35675
r52         21        -57         79        216         79   68.41144
r53          3         -6          8         23          8   7.463587
r54          1         -3          4         12          4  4.2281015
r55          1         -3          5         14          5  4.2623042
r56         14        -38         37        113         37  38.218651
r57         20        -70         48        167         48  59.329445
r58         18        -57         44        146         44  51.224711
number of observations will be reset to 58
Press any key to continue, or Break to abort
obs was 0, now 58
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[58,6]
            c1         c1         c2         c1         c2         c3
 r1        657       1092       1500       1908       1500  203.92078
 r2         48         44         65         86         65   10.32214
 r3         66         45        107        169        107  31.134122
 r4        143        153        332        511        332  89.722657
 r5        235        128        558        988        558  215.06483
 r6         43         27         66        105         66  19.448079
 r7         72         -1        201        404        201   101.4918
 r8         39         31         66        101         66  17.281394
 r9        274       -502       1385       3273       1385  944.24883
r10         40         25         58         91         58  16.451828
r11       1480       2373       2942       3511       2942  284.32859
r12        296        372       1027       1682       1027  327.62489
r13         57         32        172        312        172  69.998518
r14       3220       4876       6875       8874       6875  999.67511
r15        264        376        664        952        664  143.80225
r16         87        -54        231        517        231  143.12133
r17        117         18        230        442        230  106.06295
r18         77          8        120        232        120  55.967842
r19        363         73       1047       2021       1047  486.84331
r20        777       1005       1640       2275       1640  317.46049
r21        112         92        174        256        174  41.047275
r22        136         99        169        239        169  34.900503
r23        399        423        672        921        672  124.66943
r24        212        112        304        496        304  95.994667
r25       1410       2238       2468       2698       2468  115.20541
r26        426        509        795       1081        795  142.75727
r27        240        277        420        563        420  71.346301
r28        254      -1243       1783       4810       1783  1513.6927
r29        151        223        319        415        319  48.060843
r30        264        240        426        612        426  93.160805
r31        448        592       1049       1506       1049  228.57415
r32        597        487       1558       2629       1558  535.54705
r33        312        457        780       1103        780  161.25425
r34        184        121        277        433        277   77.93561
r35        785        927       1205       1483       1205  139.10257
r36       1124       1594       2274       2954       2274  340.07443
r37        733        929       1238       1547       1238  154.69418
r38         31         32         64         96         64  15.845961
r39         55         58         80        102         80  10.774061
r40        392        581        810       1039        810  114.52732
r41        145         95        412        729        412  158.65774
r42        156        146        234        322        234  43.770001
r43        121        116        244        372        244  64.033098
r44         51         45        115        185        115  34.758927
r45        891       1000       1527       2054       1527  263.67731
r46        204        326        485        644        485  79.276258
r47        430        576        800       1024        800   112.0124
r48        516        423        749       1075        749  163.01209
r49         36         30         88        146         88  29.237969
r50        123         79        207        335        207  64.124283
r51        572        285       1157       2029       1157  436.07774
r52        448        436       1400       2364       1400  481.79282
r53        143        140        194        248        194  26.864619
r54         33         30         51         72         51  10.417559
r55         88         80        128        176        128  24.204341
r56        171        145        288        431        288  71.468652
r57       1081       1238       3074       4910       3074  918.12087
r58        125       -286        545       1377        545  415.78217
number of observations will be reset to 58
Press any key to continue, or Break to abort
obs was 0, now 58
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. 
. 
. 
. 
. /* Report regions */
. run "$code/regions"

. foreach k in mao_mit_est mao_mit_sl mao_mit_ot mao_mit_tot{
  2. matrix mao=`k'
  3. 
. matrix mao=region'*mao
  4. scalar nregion=rowsof(mao)
  5. 
. matrix smao=mao[1...,3]
  6. matrix bmao=mao[1...,2]
  7. matrix nmao=mao[1...,1]
  8. 
. *matrix tmao=bmao
. 
. matrix cil_mao=J(nregion,1,0)
  9. matrix ciu_mao=J(nregion,1,0)
 10. forval i= 1/8{
 11. matrix smao[`i',1]=sqrt(smao[`i',1])
 12. *matrix tmao[`i',1]=bmao[`i',1]/smao[`i',1]
. matrix cil_mao[`i',1]=int(bmao[`i',1]-signi*smao[`i',1]+0.5)
 13. matrix ciu_mao[`i',1]=int(bmao[`i',1]+signi*smao[`i',1]+0.5)
 14. }
 15. 
. 
. matrix mao=nmao,cil_mao,bmao,ciu_mao,bmao,smao
 16. matrix list mao
 17. 
. svmat mao, names(v)
 18. ren v1 n
 19. ren v2 cl
 20. ren v3 e
 21. ren v4 cu
 22. ren v5 em
 23. ren v6 s
 24. 
. export excel n cl e cu em s using "$excel/${name}smi.xls", sheet("r`k'") sheetmodify firstrow(varlabels) missing(".")
 25. clear
 26. 
. 
. }

mao[8,6]
           c1         c1         c2         c1         c2         c3
c1       5667       9322      12647      15972      12647  1662.6935
c2        961       1781       3424       5067       3424  821.29419
c3       1733       3332       4313       5294       4313  490.43654
c4        929       1811       2334       2857       2334  261.63082
c5        644        936       1184       1432       1184  124.23889
c6        282       -560       1315       3191       1315   938.1697
c7        927       1936       2655       3374       2655  359.35781
c8      11143      23503      27872      32241      27872  2184.6712
number of observations will be reset to 8
Press any key to continue, or Break to abort
obs was 0, now 8
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[8,6]
           c1         c1         c2         c1         c2         c3
c1       4649       7448       8944      10440       8944  747.77845
c2       1706       1591       3084       4577       3084  746.37178
c3       1578       1102       2551       4000       2551   724.6585
c4        961       1075       1377       1679       1377  151.19449
c5        877        698       1279       1860       1279  290.73311
c6         70        -56        156        369        156  106.27689
c7        527        428        950       1472        950  260.99707
c8      10368      15637      18341      21045      18341  1352.0995
number of observations will be reset to 8
Press any key to continue, or Break to abort
obs was 0, now 8
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[8,6]
           c1         c1         c2         c1         c2         c3
c1         83         68        120        172        120   26.11062
c2        133        -44        448        941        448  246.66171
c3         91       -554        728       2011        728  641.45692
c4         20        -59         31        122         31  45.596797
c5          9        -11         20         52         20  15.796943
c6          9         -1         45         92         45  23.534714
c7         98        -22        244        511        244  133.73891
c8        443        231       1636       3041       1636  702.68011
number of observations will be reset to 8
Press any key to continue, or Break to abort
obs was 0, now 8
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

mao[8,6]
           c1         c1         c2         c1         c2         c3
c1      10399      18064      21711      25358      21711  1823.2948
c2       2800       4682       6956       9230       6956   1136.854
c3       3402       5422       7592       9762       7592  1084.9539
c4       1910       3131       3742       4353       3742    305.597
c5       1530       1850       2483       3116       2483  316.56056
c6        361       -372       1516       3405       1516  944.46337
c7       1552       2921       3849       4777       3849  463.83575
c8      21954      42522      47849      53176      47849  2663.5917
number of observations will be reset to 8
Press any key to continue, or Break to abort
obs was 0, now 8
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. 
. *************************************************************************************
. *       Differences
. *************************************************************************************
. 
. /* Report Total */
. matrix mao=sum_mit_all

. matrix dmao=mao[2,1..2]-mao[1,1..2]    /*d*/

. matrix smao=sqrt(mao[2,3]+mao[1,3])   /*t*/

. matrix tmao=dmao[1,2]/smao[1,1]    /*t*/

. matrix pmao=1-normal(abs(tmao[1,1]))    /*p*/

. matrix cil_mao=int(dmao[1,2]-signi*smao[1,1]+0.5)

. matrix ciu_mao=int(dmao[1,2]+signi*smao[1,1]+0.5)

. matrix smao=int(smao[1,1]+0.5)   /*t*/

. matrix dmao=dmao,smao,tmao,pmao,cil_mao,ciu_mao

. matrix list dmao

dmao[1,7]
            c1          c2          c1          c1          c1          c1          c1
r1        -775       -9531        2569  -3.7096666   .00010377      -14668       -4392

. 
. svmat dmao, names(v)
number of observations will be reset to 1
Press any key to continue, or Break to abort
obs was 0, now 1

. ren v1 n

. ren v2 e

. ren v3 s

. ren v4 t

. ren v5 p

. ren v6 cl

. ren v7 cu

. format %7.4f  p

. export excel n e s t p cl cu using "$excel/${name}smi.xls", sheet("dmao") sheetmodify firstrow(varlabels) missing(".")
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. clear

. 
. 
. /* Report Strata */
. matrix dmao=mao_mit_sl[1..58,1..2]-mao_mit_est[1..58,1..2]    /*d*/

. matrix smao=mao_mit_sl[1..58,3]+mao_mit_est[1..58,3]   /*v*/

. 
. matrix tmao=J(nstrata,1,0)

. matrix pmao=J(nstrata,1,0)

. matrix cil_mao=J(nstrata,1,0)

. matrix ciu_mao=J(nstrata,1,0)

. forval i= 1/58{
  2. matrix smao[`i',1]=sqrt(smao[`i',1])
  3. matrix tmao[`i',1]=dmao[`i',2]/smao[`i',1]
  4. matrix pmao[`i',1]=1-normal(abs(tmao[`i',1]))    /*p*/
  5. 
. matrix cil_mao[`i',1]=int(dmao[`i',2]-signi*smao[`i',1]+0.5)
  6. matrix ciu_mao[`i',1]=int(dmao[`i',2]+signi*smao[`i',1]+0.5)
  7. matrix smao[`i',1]=int(smao[`i',1]+0.5)   /*t*/
  8. }

. 
. matrix dmao=dmao,smao,tmao,pmao,cil_mao,ciu_mao

. matrix list dmao

dmao[58,7]
             c1          c2          c3          c1          c1          c1          c1
 r1        -243        -785         201   -3.904838   .00004714       -1186        -382
 r2         -40         -43           6   -7.095765   6.432e-13         -54         -30
 r3         -10           9          31   .29393813   .38440259         -51          70
 r4         -18         -76          89  -.85636922   .19589679        -252         101
 r5          38         -47         214  -.21963906   .41307613        -474         381
 r6          -5          10          19   .51482976     .303336         -28          49
 r7         -31         -95          99  -.95938199   .16868317        -292         103
 r8         -18         -10          11  -.94148693   .17322769         -30          11
 r9        -154       -1106         944  -1.1714974   .12069946       -2993         782
r10           2          18          16   1.0941034   .13695482         -14          51
r11        -136       -1293         281  -4.5961607   2.152e-06       -1855        -729
r12        -133        -739         303  -2.4392504   .00735888       -1344        -132
r13          -2         -39          50  -.78502991   .21621801        -137          60
r14         -72       -1270         770  -1.6492964   .04954347       -2809         270
r15          45         -26         141  -.18419084   .42693187        -307         256
r16          26         -73         143  -.51017507   .30496441        -358         213
r17           7          -3         106  -.02829144   .48871485        -214         209
r18          35          54          56   .96520239   .16722173         -57         166
r19          -5        -521         487  -1.0701619   .14227321       -1494         453
r20           1        -445         317  -1.4018131   .08048554       -1079         190
r21         -39         -54          41  -1.3172415   .09387884        -135          28
r22          -5          22          35   .63084876   .26406971         -47          92
r23         -28        -132         125    -1.05889   .14482493        -380         117
r24          87         161          96   1.6775932   .04671327         -30         353
r25        -541        -472         115  -4.1014467   .00002053        -701        -241
r26          55         105         143   .73582562   .23091841        -179         390
r27          56         174          71   2.4391766   .00736039          31         317
r28         -71       -1478        1514  -.97642367   .16442728       -4504        1549
r29          51         155          48   3.2251852   .00062946          59         251
r30          60         153          93   1.6434994   .05013982         -32         339
r31          92        -121         229  -.52942892   .29825397        -577         336
r32        -221         157         535   .29320177   .38468397        -913        1228
r33         103         331         161   2.0612516   .01963952          10         652
r34          53          -6          78  -.07703092   .46929947        -161         150
r35        -205        -471         139  -3.3950392   .00034309        -747        -193
r36        -194        -857         340  -2.5201027   .00586603       -1536        -176
r37        -319        -516         154  -3.3426394   .00041493        -824        -206
r38          -3          -4          16  -.25255728   .40030517         -35          28
r39         -12          11          11   1.0220718   .15337345         -10          33
r40          89         416         114    3.634316   .00013936         187         645
r41         -76        -205         159  -1.2923576   .09811666        -521         112
r42          32          86          44   1.9652761   .02469115          -1         174
r43         -49        -101          64  -1.5780999   .05727133        -228          27
r44         -37         -75          35  -2.1579513   .01546581        -144          -4
r45          19        -202         263   -.7673377   .22144038        -727         324
r46          56          47          79   .59290375   .27662279        -111         206
r47         168         336         110   3.0459325    .0011598         115         557
r48         165         214         163   1.3133532   .09453199        -111         540
r49         -16         -41          29  -1.4047342   .08005017         -98          17
r50          49          83          64   1.2944102   .09776188         -44         211
r51         192         510         374   1.3638712   .08630406        -237        1258
r52         -15        -517         477  -1.0840595   .13916921       -1470         437
r53         -70         -66          26   -2.557443   .00527224        -117         -13
r54         -12         -13          10  -1.3654085   .08606235         -31           6
r55         -21          -7          24   -.2937955   .38445709         -54          41
r56           7          25          60   .41396727   .33944904         -95         146
r57         619        -246         916  -.26849977   .39415733       -2077        1586
r58         -81        -453         413  -1.0978767   .13612917       -1277         372

. 
. svmat dmao, names(v)
number of observations will be reset to 58
Press any key to continue, or Break to abort
obs was 0, now 58

. ren v1 n

. ren v2 e

. ren v3 s

. ren v4 t

. ren v5 p

. ren v6 cl

. ren v7 cu

. format %7.4f  p

. export excel n e s t p cl cu using "$excel/${name}smi.xls", sheet("sdmao") sheetmodify firstrow(varlabels) missing(".")
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. clear

. 
. 
. 
. 
. /* Report regions */
. matrix mao_mit_sl=region'*mao_mit_sl

. matrix mao_mit_est=region'*mao_mit_est

. 
. matrix dmao=mao_mit_sl[1..8,1..2]-mao_mit_est[1..8,1..2]    /*d*/

. matrix smao=mao_mit_sl[1..8,3]+mao_mit_est[1..8,3]   /*v*/

. 
. matrix tmao=J(nregion,1,0)

. matrix pmao=J(nregion,1,0)

. matrix cil_mao=J(nregion,1,0)

. matrix ciu_mao=J(nregion,1,0)

. forval i= 1/8{
  2. matrix smao[`i',1]=sqrt(smao[`i',1])
  3. matrix tmao[`i',1]=dmao[`i',2]/smao[`i',1]
  4. matrix pmao[`i',1]=1-normal(abs(tmao[`i',1]))    /*p*/
  5. 
. matrix cil_mao[`i',1]=int(dmao[`i',2]-signi*smao[`i',1]+0.5)
  6. matrix ciu_mao[`i',1]=int(dmao[`i',2]+signi*smao[`i',1]+0.5)
  7. matrix smao[`i',1]=int(smao[`i',1]+0.5)   /*t*/
  8. }

. 
. matrix dmao=dmao,smao,tmao,pmao,cil_mao,ciu_mao

. matrix list dmao

dmao[8,7]
            c1          c2          c3          c1          c1          c1          c1
c1       -1018       -3703        1823   -2.031147   .02112004       -7348         -56
c2         745        -340        1110   -.3063691   .37966182       -2559        1880
c3        -155       -1762         875   -2.013671   .02202204       -3511         -11
c4          32        -957         302  -3.1670263   .00077003       -1560        -352
c5         233          95         316   .30047491   .38190747        -536         727
c6        -212       -1159         944  -1.2275331   .10981115       -3046         729
c7        -400       -1705         444  -3.8389072   .00006179       -2592        -816
c8        -775       -9531        2569  -3.7096666   .00010377      -14668       -4392

. 
. svmat dmao, names(v)
number of observations will be reset to 8
Press any key to continue, or Break to abort
obs was 0, now 8

. ren v1 n

. ren v2 e

. ren v3 s

. ren v4 t

. ren v5 p

. ren v6 cl

. ren v7 cu

. 
. format %7.4f  p

. export excel n e s t p cl cu using "$excel/${name}smi.xls", sheet("rdmao") sheetmodify firstrow(varlabels) missing(".")
file C:/Users/silvio/Documents/CVR/ryp/excel/s2_dokrigsmi.xls saved

. 
. clear

. capture log close
